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Introduction 



Study of spectral fluctuations in a quantum dynamical system was initi- 
ated by Wigner in 1956 [TJ . He provided explanation of observed shortage of 
close spacing in nuclear energy levels using statistical arguments. For exam- 
ple, one can obtain precise information about levels ranging from number N 
(of order 10^) to + n from neutron-capture experiments in heavy nuclei. 
Theoretically, it is difficult to analyse levels of this order on the basis of shell 
structure and collective or individual particle quantum numbers. Wigner 
successfully shown that the highly exited states can be understood by as- 
suming a working hypothesis that all shell structure is washed out and that 
no quantum numbers other than spin and parity remains good via statistical 
theory of energy levels. The aim of statistical theory is not to predict the 
detailed sequence of levels, but to describe the general appearance and the 
degree of irregularity of the level structure in the given quantum system. 

This view led Wigner to surmise a possible spacing distribution based 
on the assumption that matrix elements of the Hamiltonian matrix were 
unknown and unknowable and are randomly distributed. Success of Wigner 
surmise in explaining the spacing between energy eigenvalues resulted in 
further development of a one of the tools to study statistical properties of 
levels, known as Random Matrix Theory (see 

The random matrix theory of Hamiltonian systems is based on the as- 
sumption that we know very little about the dynamics of the given system 
except certain symmetry properties. These symmetry properties impose re- 
strictions on the form of the Hamiltonian matrix. Based on the symmetry 
properties a normalised probability distribution on the elements of the Hamil- 
tonian matrix can be obtained using a metric in the space of matrix elements 
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which is invariant under a similarity transformation. The joint probabihty 
distribution is chosen to minimize information about the Hamiltonian matrix 
and subjected to the condition that it is normahzed to one and has elements 
that remain finite. 

For real symmetric Hamiltonians (e.g. systems having time reversal sym- 
metry,), the similarity transformation is orthogonal and the probability dis- 
tribution is said to be described by a Gaussian Orthogonal Ensemble (GOE) 
of Hamiltonian matrices. For other interesting cases such as Hermitian and 
real quaternion Hamiltonians, the similarity transformations are unitary and 
symplectic respectively. The probability distribution of these is said to be 
described by Gaussian unitary ensemble (GUE) and Gaussian Symplectic 
Ensembles (GSE). 

There are number of statistical properties of eigenvalues of random ma- 
trices that are commonly used in the analysis of spectral properties. For, ex- 
ample spacing distribution (or nearest neighbour spacing), the A3— statistic 
(or spectral rigidity) of Dyson and Mehta [21 E] ■ 

Random matrix theory (RMT) have been successful in explaining univer- 
sality in spectral fluctuations (up to second order correlation) for different 
classes of chaotic systems as above in an appropriate high energy limit. RMT 
is expected to be applicable only on those time scales where the variables as- 
sociated with the classical dynamics are random enough, to fully randomize 
the matrices, associated with the corresponding quantum operators. More- 
over, due to underlying assumptions on which framework of RMT is build, 
it is not expected to shed any light on non-universal behaviour of spectral 
fluctuations which is the characteristic of a given system. 

A new dimension is added to the study of spectral properties of a quantum 
system after the conjecture of Percival j3] that there exist a regular and an 
irregular quantum spectra corresponding to the integrable an non integrable 
classical dynamics. This indication that the nature of classical dynamics 
can have its imprints in the behaviour of a corresponding quantum system 
provided a new direction to study of statistical properties of levels. 

Classical Hamiltonian systems can display wide variety of dynamical fea- 
tures from integrability to chaos. Before the works of Poincare, classical dy- 
namics was viewed as a paradigm of regularity and predictability. The works 
by Poincare, and subsequently by Lyapunov, BirkhofF and others opened the 
complicated world of classical dynamics. Within a few years after these con- 
cept of determinism was abandoned from classical mechanics. The concept 
about classical motion changed radically and led to birth of what is nowadays 



4 



called chaos. Broadly now, classical systems can be classified as one showing 
complete integrability and others showing either complete chaos or mixed 
behaviour. 

For integrable (or regular) systems all constants of motion exist every- 
where in the phase space. A typical trajectory resides on a torus like constant 
energy surface. On other hand one in case of completely chaotic systems no 
constant of motion exist and a typical trajectory will occupy whole phase 
space available. Most of the classical systems, however, belong to an interme- 
diate regime between these extremes, showing a mixed behaviour. Typically, 
some of the trajectories will show regular behaviour winding around tori. 
These trajectories are separated by the other orbits that explore constant 
energy surface ergodically. 

There are also many Hamiltonian systems that are just step away from the 
complete integrability, known as pseudo-integrable systems. For these sys- 
tems constants of motion do not exist on a countable set of singular points. 
These systems have /-dimensional (/ being dimension of the system) sphere 
with g handles as an invariant surface embedded in 2/-dimensional phase 
space. Typical examples of pseudo-integrable systems are rational of polygo- 
nal billiards, where a particle moves freely inside polygonal enclosure whose 
each angle is of form ir/n, reflecting specularly from the walls. Results of 
some recent studies indicate that generally one can expect coexistence of al- 
most integrable and almost chaotic regions in the mixed phase space. In this 
work we are concerned with spectral properties of these systems. 

The imprints of classical behaviour on a quantum system can best be un- 
derstood via semiclassical techniques which provides a information regarding 
quantum dynamics in terms of properties of corresponding classical system. 
Semiclassical techniques had been developed since the beginning of the quan- 
tum mechanics. Due to easiness with classical mechanics the search for con- 
nections between classical and quantum mechanics has been naturally sought 
for. The WKB approximation was the first method built to extract quantum- 
mechanical properties using information of classical orbits. The recent inter- 
est in non-integrable systems further initiated development of semiclassical 
techniques which are applicable for large class of classical systems. 

In this work we will consider dynamical systems as mentioned above 
known as pseudointegrable billiards. Due to mathematical intractability of 
these systems, few exact results are known. Numerical studies performed on 
the energy spectra of these systems led to the belief that the spectral statis- 
tic is intermediate between those for the Poisson ensemble and the GOE of 
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random matrices. However, the results on various measures of the spectral 
statistics as analysed from the periodic orbit theory have been indicative only 
and do not bring out a complete or explicit picture of underlying correlations. 
Moreover, numerical studies have to deal with finite number of levels only. 
There is thus possibility that numerical results may not represent asymptotic 
trend. 

Our aim here is to extend semiclassical formalism to obtain explicit ana- 
lytical expressions for various important measures of spectral statistics. This 
will therefore can bring out dependence of spectral properties on different 
system dependent parameters. 



Plan of the thesis is as follows: 

In chapter one we will discuss main issues involved in development of 
proper semiclassical framework that can be applicable for wide class of sys- 
tems. The problem with WKB or similar methods developed earlier is that 
they were applicable only for systems that are separable. Existence caustics 
prohibits these methods to be useful beyond short time scale. Major im- 
provement in semiclassical methods has been taken place via Path Integrals 
which ultimately leads to a semiclassical approximation to density of states 
in terms of summation over periodic orbits of the classical system. 

As the periodic orbits and their attributes such as stabihty, action etc. 
becomes important in establishing a semiclassical approximation, study of 
the nature of classical system in question becomes mandatory. In the case 
of polygonal billiards there exist a large body of knowledge about classical 
behaviour in the literature. Some of these are very helpful in the study of 
periodic orbits and their various properties that are relevant to semiclassical 
methods. These will be discussed in chapter two. 

In chapter three we take up the first task [HI: detailed study of periodic 
orbits in some pseudointegrable billiards. Exploiting geometrical properties 
of the billiards considered we develop a modified version of interval exchange 
maps, called as polar maps. From these maps we show how to classify and 
enumerate all the periodic orbits of these billiards. We will also quantify 
actions and other parameters required for the semicalssical study. 

In chapter four a very important attribute of periodic orbit i.e. growth 



6 



rate (with actions, periods or lengths) of periodic orbits is studied 0. We 
obtain exact (asymptotic) law for some pseudointegrable billiards. We show 
that this law is quadratic in length and determine constants associated with 
it. Some generalization of this law is also suggested. 

Using information gathered in chapter 3 and 4 we obtain explicit analyt- 
ical expression for two point correlation function and its fourier transform in 
chapter five. Generally semicalssical study ends after establishing a sum over 
periodic orbits. This does not bring out dependence of spectral fluctuations 
on the various attributes of periodic orbits explicitly. We here, go further 
to carry out this summation by converting it into a integration via proper 
measure, which enables us to obtain analytical forms for spectral fluctuation 
measures. This also brings out role played by growth rate of periodic orbits 
explicitly. 

In chapter six we use two point correlation function to obtain short and 
intermediate range spectral measures such as spacing distribution, number 
variance and spectral rigidity [7j. We will also compare our results with some 
numerical experiments. 

In the concluding chapter we discuss our results and outline some of the 
future research activities in this field. 
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Chapter 1 



Modern Approaches to 
Semi-Classical Quantization 

1.1 Introduction: 

Semi-classical methods have been developed since the beginning of quantum 
mechanics because of the importance of the models based on classical par- 
ticles and fields. Bohr and Sommerfeld initiated these methods to quantize 
atomic systems. However their methods implicitly assumed separability be- 
tween the different internal degrees of freedom and hence were of limited 
utility since most of the classical systems are not separable[T]. For exam- 
ple even the three-body helium atom problem can not conform to the Bohr- 
Sommerfeld scheme without further assumptions. There were many attempts 
to resolve this problem e.g. see [2|. However, major development in this di- 
rection had taken place with the advent of Heisenberg's matrix mechanics 
and of Schrodinger's wave mechanics. 

The Schrodinger wave equation of quantum mechanics turns out to be 
extremely useful in the development of semi-classical theory since it allows 
a discussion of the relation to classical mechanics along the lines similar to 
the relation between wave optics and geometric optics, albeit in the multidi- 
mensional phase space rather then in the three dimensional physical space. 
This similarity was developed by Brillouin, Wentzel and Kramers for one 
degree of freedom (well known as WKB approximation) based of well devel- 
oped mathematical framework of asymptotic methods of linear mathematical 
physics 13 El . For larger systems that are separable and reducible to several 
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one-degree of freedom systems variant of this method (EWKB approxima- 
tion) was developed. 

Since semi-classical methods carry out the quantization on the basis 
of knowledge of the classical motion, they are susceptible to the difficul- 
ties caused by the non-linearities of Hamilton's equations of classical me- 
chanics. The inapplicability of standard WKB methods to non-linear sys- 
tems prompted further development of semi-classical methods |12[ ll3j. These 
methods known as the periodic orbit quantization have been applied to a 
large variety of systems. These methods basically use a semi-classical ap- 
proximations of the trace of the resolvent of the quantum Hamilton operator. 
This trace formula, which establishes link between the quantum operator and 
classical periodic orbits, can be used to obtain approximate values of quan- 
tum eigenenergies for wide class of systems. The trace formula which contains 
different terms depending on the nature and on the stability of the periodic 
orbits. It not only improves the accuracy of semi-classical methods but also 
provides a framework for the description of diffraction, non-linear stability 
effects, bifurcation of periodic orbits etc. 

In this thesis we will discuss the imprints of nature of classical Hamilto- 
nian system on the corresponding quantum mechanical system, in particular 
on the collective properties of quantal eigenenergies. We will use the frame- 
work of semi-classical methods for our study. Before coming to the main 
theme, we will discuss in this chapter some important issues in modern semi- 
classical methods. 



1.2 Short- wavelength Asymptotics of Schrodinger 
Equation: 

Consider the Cauchy problem with rapidly oscillating initial data for the 
Schrodinger equation [U [7j 

V'|i=o=^(q)exp(^5o(q)) (1.2) 

where q G -R", functions V^(q), >S'o(q) are real valued and infinitely differen- 
tiable, Ao(q) is infinitely differentiable with compact support. One can seek 
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an asymptotic solution of the problem as ^ — Oand q G -R", < t < T. 
The equation (jl.2p can be regarded as the leading term of an asymptotic 
expansion in h, in which 0(^)terms are neglected. It is also assumed that 
an asymptotic solution of equation (jLlj) 

^(q,t) =A(q,t)expQ^(q,t)) (1.3) 

which is valid only for the short elapsed time t, after which it should be 
replaced by the sum of higher order terms in an asymptotic expansion. The 
reason for this single term break down is formation of caustic which occur 
at time of order (h^). However for sufficiently longer time above assump- 
tion breaks down altogether. We will elaborate these issues in the following 
paragraphs. 

Substituting equation ()1.3|) in the Schrodinger equation and neglecting 
terms of order (h) and higher, we get the time-dependent Hamilton- Jacobi(H- 
J) equation for the action S (q, t), 

where H is Hamiltonian and the canonical momentum p of iJ(q, p,t) is 
replaced by p = dS/dq. The order (h) term gives amplitude transport 
equation for A or the continuity equation for p(q, t) = |A(q, )f:)p as 

^^ + ^[pi^,t)v{q,t)] = (1.5) 

where the velocity field v = dH{ci,p,t)/dp with p = dS/dq. One then 
solves equation ()1.4|) for S{q,t) and use this S to solve equation ()1.5|) . 

The H-J equation is typical of non-linear equations, hence has a bewilder- 
ing variety of solutions. Among them are the complete integrals which gives 
all the trajectories of the system. As we are concerned with those solutions 
that have a relationship to a quantum mechanical wave function, it is desir- 
able to find a way of defining S such that it will be a single valued function 
of its variables. The first step to do this is to use phase space description of 
trajectoriesjSl E]- The concept of Lagrangian Manifold becomes important 
here. 
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1.2.1 Lagrangian Manifold: 

The initial action 5'o(q) and momentum field po = dSo/dq can be viewed as 
a vector field on the /-dimensional configuration space. This field represents 
the initial momenta of the swarm of particles whose initial density is given 
by Po(q) = |^o(q)P- This swarm of particles and vector field constitute 
the classical and semi-classical interpretation of the initial action function 
in configuration space. Since the momentum field imposes / independent 
constraints on the 2/ variables (q, p) of the phase space, initial swarm of 
the particles lies on the / -dimensional surface Aq. The surface Aq is graph 
of the function p = po(q) i.e. the set of points in phase space of the form 
(q, Po(q)). This surface satisfy properties of Lagrangian manifold, A, the 
formal definition of which is as follows: 

Definition 1 Lagrangian manifold is a f-dimensional surface A in 2f- dimensional 
phase space such that V(q, p) G A and for W 6z = (5q, 6p) tangent to A at 
(q, p) representing a small displacement , the action of the symplectic form 
(i.e. differential 2- form) on all 6zi,Sz2 G {Sz} defined as 

w{5zi, 5z2) = 6zi ■ JT^ ■ 5z2 = 5pi • 5q2 - 5p2, 5c[i = (1.6) 

where J_ is unit symplectic matrix. 

All the curves in the 2-dim. phase space are Lagrangian manifolds. The 
concept is therefore really needed for multidimensional problems. It can be 
immediately seen that surfaces q = const, or p = const, are A. Since the 
value of symplectic form of ()1.6p is invariant under canonical transformations 
every A is a constant q or constant p surface in some set of canonical co- 
ordinates. The graph of any curl- free momentum field is surface A. The 
converse is however, not true, since A may contain points at which derivatives 
in the curl-free condition {i.e.d'Pi/ d(\j) are not defined. Such points are 
generally associated with caustics. 

The /-dimensional vectors q and p can be regarded as smooth functions 

of M = (mi M2, j""/); " co-ordiuates" labelling of A. The variables c\{u) are 

locally invertible if Jacobian |5q/(9p| do not vanish. Then and only then one 
can define p(q) = p(-u(q)), a function of q on A. If above Jacobian vanishes 
then 

dpi_ dpi_ duk 
dqj ^ duk dqj 
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will behave badly. In 1— dimension since dpi/du^ 7^ it will diverge. In 
higher dimensions some of the eigenvalues of [dp/du] may vanish at the 
same place where some of the eigenvalues of [9q/ du] also vanish, resulting 
in the complicated behaviour. A phase space vector 6z tangent to A can be 
written in the u— co-ordinates as dz = (^^6u, ^6uj. If the matrix [dq^/du] 
is singular, then for all 6u ^ and 6q = , 6z has vanishing components 
in the g-components resulting in a caustic in the configuration space. The 
order of caustic or the number of null eigenvectors 6u, is given by co-rank of 
the matrix [dq/du]. The order of caustic as well as set of points for which 
dq/du = are independent of the choice of u on A. The important theorem 
about A is 

Theorem 1 A region of A which is free of singular points if projected on 
configuration space gives curl-free momentum field. ^ 

Since flow preserves the symplectic structure (cf. Liouville Theorem) 
A evolves into another A under time evolution. Any function 5'(q) which 
satisfy p(q) = dS/dci on A will now be called generating function of A. 

Theorem 2 On a region of A, free of singular points, one can define unique 
generating function up to additive constant.^ 

This additive constant is usually associated with phase conventions in 
semi-classical applications. For A, having region of singular points, one can 
divide A into sub-regions which extend up to and separated by the caustics 
as shown in Fig. 1.1. Each such region corresponds to a distinct branch of the 
momentum fleld p;,(q) and associated generating function S'b(q) up to its own 
additive constant. One can deflne action function ^(q, p) on A itself, as the 
line integral J p ■ dq along the contour belonging to the A and taken relative 
to an arbitrary initial point, see Fig. 1.2. This can be done by demanding 
that the different functions 5*5 (q) approach one another at the caustic di- 
viding branches. In this way one can also link some or all of the additive 
constants together. The function ^(q, p) will be multi- valued not only due 
to caustics but also due to non-trivial (e.g. multi-connectedness) topology of 
A. The generating function in general is speciflc to the choice of co-ordinates 
Some A consist entirely of caustic points (e.g. x = const, in 1-dimension). 
Such manifolds do not have generating function S{x). However, with respect 
to other co-ordinate systems one can define a generating function for such 
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Fig.1 .1 : The existance of a caustic point {\,Pq) on 
a Lagrangian manifold A, leads to a multivalued 
momentum field, p=p(x), shown here by points 1,2 
for one deg ree of freedom. 
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Fig. 1.2: An action function S(x,p) can be defined at the 
pfiase space points ona Lagrangian manifold along 
a contour confined to the manifold, relative to point "0". 
This action function becomes multivalued upon projection 
onto configuration space, but is continuous at caustics. 
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surfaces. The locations of the caustic points are relative to the representa- 
tion being used. This fact plays important role in the WKB theory. It is 
possible to cover every A (Fig. 1.3) with overlapping regions such that every 
region is caustic free in some representation obtained from a commutating 
mixture of q's and p's. A whole family of wave functions, such as a complete 
set of commutating observables, will reproduce a whole family of A mani- 
folds parameterized by some set of parameters A = (Ai^A2,A3^ Xf). Thus 
phase space is foliated up into /- parameter family of /-dimensional A man- 
ifolds. The function S{q, A) is one of the generating functions of canonical 
transformations in classical mechanics. 

1.2.2 Solution of H- J and Amplitude Transport Equa- 
tions: 

Now consider an initial Lagrangian manifold A', free of singular points, ob- 
tained from initial action 5'(q',t') = 'S'o(q') and p(q',t') = dS{c{ ,t')/dc{ . 
Each point of A evolves under Hamilton's equations, mapping into another 
A at time t . If A is also free of singular points then one can define a 
generating function S{q",t"), which is also a solution of H-J equation, as 

S{c^',t")^S{c^,t') + R{c^',t":c^,t') (1.7) 

q ,t 

where -R(q", t" : q', t') — J [p ■ dq — Hdt] is Hamilton's principal function. 

q ,t 

As t" — )• t\one gets R — and S{q^',t") — S{q',t') thus satisfying 
required initial conditions. 

Since the particle density is conserved one can write p(q", t")dc[ — p(q', t')dci , 
and since p = |y4p, we get 

Ai^',t")=A{^,t') 

Here the absolute sign in this equation required only when one considers 
many branches. 

Thus the WKB solution to the Cauchy problem can be written as 

*(q",O = ^(q',0 

(1.9) 



det 



(9q 



9q" 



(1.8) 



det 



dq' 



exp{l[5(q',0 + W,^":q',0] 
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X 



Fig. 1.3: In one degree of freedom, configuration space 
caustics (2&2') and momentum space causics (l&V) never 
occur at the same place. In higher degrees of freedom also 
there always exist a representation in which a given point is 
caustic free. Thus a Lagrangian manifold can always he 
covered hy overlapping regions (e.g. A,B,C,D) that are 
caustic-free. 
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This expression diverges when q" is at caustic (i.e. det9q'/9q" diverges) of 
A". This is illustrated in Fig. 1.4 for one degree of freedom. This divergence 
represents a non uniformity in the variables (q", t") of \E' in since the error 
in ^, — > 0, like fi, as Ti — > for fixed (q",t"), however it goes to infinity 
for fixed h as (q",t") approaches caustic. Though the solution near caustic 
is not valid one can continue the solution through the divergence. However, 
function q'(q",f,f') is then multi- valued. By passing through caustic the 

determinant changes sign via divergence. Thus det becomes imagi- 

nary after first caustic pass. One can absorb this imaginary part into a phase 
factor as exp(— i?77r/2), forcing A to be positive. The 7] here is well known 
as Maslov index. Considering multiplicity of orbits arriving at q" at time t'\ 
one can write 

*(q", t") = E A(q", t") exp {^5,(q", t") - zr7,|} . (1.10) 

To compute indices r]'s one needs to consider momentum-space wave func- 
tion which is Fourier transform of a configuration-space function. The Fourier 
transform integral is evaluated by stationary phase approximation. This leads 
to the momentum-space function 

$(p) =exp(za0 |detM(q)r^exp|^[5(q)-q-p]| (1.11) 

where the symmetric matrix M(q) = dp/dq — d'^S{q)/dqdq and a is an 
integer given by the index of inertia of M_ (i.e. number of positive eigenvahics 
minus number of negative eigenvalues). If A do not have singular points in 
say, X — representation but do have a caustic in P — representation, then at 
the momentum-space caustic det M(q) — > 0, hence a changes discontinuously 
across the caustic by ±2 • n, since number of eigenvalues of M which vanish 
at caustic is same as the order "n" of the caustic. Therefore relative phase 
shift between the two branches will have form exp(— ir^7r/2), and 



$(p) = exp [iao^ ^ A (p) exp |^ 



;i.i2) 



where ckq is index of inertia of one of the branches for which = 0, and 
Sb{p) = 5'6(q) - q • p, Afe = (p6(p))^''^ with ^^(p) = p6(q)| det 9q/9p|. 

Thus the relative phase shifts between branches can be determined by 
switching to a representation which is caustic free. Considering any two 
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X 



Fig. 1.4: The orbit a, has reached a caustic, b has not yet 
while c, has already passed one ■ After caustic deveoped 

X andp at corres ponding point (e.g. a') becomes 
multivalued . 
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branches the phases of the two branches of $b(p",t") can be written as?7 
Sb{p) — irib'n'/2 + iabir/'i. Demanding continuity of phases across caustic in 
caustic free representation, we get 



The bounds on Ar/, change in rj across a caustic is then —n < At] < +n, n 
being order of caustic. By extending process, it is possible to associate total 
Ai] measured between the end points of a curve segment on A, which may 
cross a number of caustic and straddle several branches. The quantity Arj 
associated with a directed curve segment of A is the Maslov index of that 
segment. The index depends only on the geometry of A and its' projection 
on configuration space. It also depends on the representation chosen. 

The method discussed in this section is simple and allows us to derive 
several important expressions. However, caustics are causing such difficulties 
that formulation has restricted applicability. The determinant det(9p/9q) 
plays an important role in the formulation which becomes singular at caustic. 
At caustic, the method of matching mentioned above uses different equations 
than preceding ones. This method is also difficult to apply for higher order 
terms in h expansion. For this reason it is advantageous to reformulate the 
problem in terms of Path Integrals. 



1.3 Path Integrals and Semi-classical Limit 



An alternate semi-classical quantization technique via Feynman Path Inte- 
gral gained momentum due to problems discussed in the previous section. 
The first hint regarding this comes from the work of SelbergjS] who obtained 
a path integral formulation for a particle moving on a Riemann surface of 
negative curvature. On the short time scales path integrals had been used to 
study semi-classical integrable svstemsP|llUj. After overcoming the problem 
with caustics the semi-classical path integral was extended to long times [TT]. 
These efforts culminated in the trace formulae of Gutzwillar, Balian and 
Bloch, as well as Berry and Tabor [T^ ITi] . The trace formulas are among 
the very few theoretical results of any generality that connects quantum me- 
chanics with classical mechanics via semi-classical approximation for quantal 
density of states of a bound quantum system, in terms of a sum over the 
periodic orbits of corresponding classical system. In this section we briefly 
outline these developments. 



V2 = Vl + 



a2 — «i 




2 
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1.3.1 The resolvent operator and its' singularities 

Consider J3] time evolution of equation which is linear in nature. 

Therefore we can define a linear operator U {t, t') such that 

vI/i(q) = f/(t,t')^Hq) (1-14) 

The operator U known as evolution operator, obeys the Schodinger equation 
p.ll) . with the initial condition U{t',t') = I, an identity operator on the 
space of state vectors. The Hamiltonian defined as 

^. = -E^V? + nq) (1-15) 
j=i ■^"h 

is a Hermitian operator {H = H'^), which implies 9t (^'t (q) l^'t (q)) = 0. 
Hence the evolution operator U is unitary (u^'U = . The formal solution 

for U{t,t') can be written as 

. i 

U{t,t') = i -I J dsHsU{s,t') (1.16) 
t' 

In case of time-independent Ht one can solve this equation by iterations to 
obtain the Dyson series for the operator which finally reduces to 

- . r. i iH(t-t')] 

U{t,t')=expl V^j • ^^-^^^ 

The evolution operator also obeys group property JJ(t,t")U(t" ,t') = (j{t^t'). 
(Using this we can set t' = 0). One can introduce evolution operators for 
' 'forward' ' and ' 'backward' ' propagation in time such as 

±e{±t)u{t) = u^{t), (1.18) 

where Q{t) is the Heaviside step function (i.e. 6(t) = IVt > 0, 6(t) = OVt < 
0). The form of (|1.17j) is indication of conjugacy of energy and time variables 
as involved in Fourier transform. Thus one can describe a system either in 
the time-domain or in the energy domain. Both of these description are 
related by a Fourier transform. 
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The Fourier transform of the U^{t) is 

r - ^ (1-19) 

where E is the eigenvalue of the stationary Schrodinger equation ff|\E'(i?)) = 
El'^lE)), obtained from taking Fourier transform of (jl.lj) . The Integral 
fll.l9|l does not exist, hence Green operator does not have meaning as an 
operator-valued functions, but can be meaningful as operator-valued distri- 
bution. We can redefine G±{E) by replacing E E ± irj,!] > (i.e. by 
analytical continuation in complex E surface), then 

G±{E)= \im{E±ir]- H)-^ (1.20) 

The object G{z) = [z — H)^^ is known as resolvent of H. The operators 
G±{E) are now called as forward and backward Green operators. The resol- 
vent in bounded, except for the values of z which corresponds to the spectrum 
of Hamiltonian. Since, eigenvalues of the Hamiltonian, being a Hermitian op- 
erator , are real and G(t) -|- Q{—t) = 1, one can write evolution operator via 
inverse Fourier transform of G{z) as a contour integral 

For t > {t < 0), the contribution of the contour G^ (C* ) is zero. The resol- 
vent has two types of singularities on the real axes. The discrete spectrum of 
H results in poles on the real axis below an energy threshold Eq. The Con- 
tinuous spectrum on other hand results in the branch points (in general of 
order two, since energy is related to momentum by relation E = Ec +p'^/2m) 
at real energies {Ec} . Thus spectral decomposition of the evolution opera- 
tor can be carried out in terms of a sum over bound states |\E'fe) and sum 
over continuous states |\E'c'(i?)), of different channels G, with eigenvalues 
exp {-^} , exp {-^} respectively, as[T7| 

OO 

U{t) = Y: exp |-^} {^,\ l dEl^i'ciE)) exp {-^} {^cm 

(1.22) 
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The eigenfunctions of If are also the eigenfunctions of H. The energy 
eigenvalues are then poles of the Green operator, and hence also of its trace 
|16j . Considering bounded systems only we can write the expression for the 
level density (or density of states) d{E) as 

oo 

d{E) ^ES{E- E^) = Tr6{E - H) 

n=Q . [i.Zo) 

— - lim„_»o , Im Tr — 

The level density is derivative of the staircase function defined as the number 
of eigenvalues below energy E : N[E) = Number {En < E} = Q{E — 
En). The singularities of the level density on the real axis give eigenvalues 
of the discrete spectrum of H. It is difficult to establish exact quantum me- 
chanical expression for the evolution operator. The path integral formulation 
is proved to be quite useful in establishing an expression for the evolution 
operator and a propagator. 

1.3.2 A Propagator and Path Integral 

In position representation equation ()1.14|) . can be written as 

*t(q) = {q\U{t,t')\%,) (1.24) 

= Jdq' {q\Uit,t')\q'){q'\^t') 
= !dq'K{qq' ■tt')^t'{q') 

where K{qq' : tt') = {q\U{t,t')\q') is propagator, which gives wave function 
at time t and position q once the wave function at t', q' is known. It may be 
noted that equation ()1.24|) is in fact the expression of Huygens principle of 
the wave propagation. The propagator satisfies Schrodinger equation 
along with initial condition 

lim K{qq' : tt') = 6{q - q). (1.25) 

It also follows from the group composition properties of the evolution oper- 
ator and definition of the propagator that the propagator also satisfies the 
semi-group property e.g. 

K{qq' : tt') = J dq"K{qq" : tt'')K{q"q' : t"t') (1.26) 
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Thus transition from {q',t') — > {Qit) involves all the possible points, hence 
all possible paths. This notion of propagation over all possible path is very 
important in the path integral formulation of quantum mechanics [18j. 

Consider a transition between two space-time points (g^, tj) and {qf,tf). 
Divide the time interval r = tf — ti m N equal pieces of length At = t/N. 
The semi-group property implies (with g^v, = tf and q^, to = qi, ti) 

K{qfqi : tfti) = ■ dqidq2 ■ ■ ■ dq^^i \{ K{qn+iqn ■ At) (1.27) 

•' n=0 

The paths entering in the above equation are not necessarily classical tra- 
jectories, and in general continuous but non-differentiable. Formally, these 
paths are similar to those of Brownian motion, which are Markovian Pro- 
cesses. Using equation ()1.17|) . the propagator for the Hamiltonian of form 
(ll.lSp . over a small time interval At can be written as 

K{qn+iqn ■■ At) = exp{-^^^}|g„). (1.28) 

n 

Using identity exp{A+B} = exp{y4} exp{5} exp | — -B]|---, where com- 
mutator [A, B] = AB — BA and dots represents higher order commutators, 
one can write 

exp{-_}^exp{-— }exp{-— }^l + -(-) [T V]j (1.29) 

In the limit — > oo. At becomes small enough so that O(At^) terms can 
be neglected. Therefore, equation p.28|) becomes 

K{qn+iqn ■■ At) = (g„+i|exp| ^||gn)exp| 



which can be rearranged in form 

^ ' ' ' ' ' m /q„+i - Qr, 



m \~2 \ iAt 



Kiq^^.qn : At) = ( j exp ^ 



2 V At 



V 



Thus in limit — > oo, equation ()1.27p becomes 



ir(q.qo : At) = Jm^ / " " " / (^) ' U % (2vraAt)-^ exp U / dt L(q, q 



u 

1.30) 
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where, L(q, q) = T — V, is the classical Lagrangian. In the phase of above 
equation we have Hamilton's Principal Function 

if 

R[q{t)] = J L(q,q)c?t (1.31) 

The physical path (or classical trajectory) will be realized only if one extr- 
mizes -R[q]. 

We are now interested in the semi-classical limit of equation ()1.30|) . The 
semi-classical limit is obtained via method of stationary phase approxima- 
tion, which amounts to an expansion about the classical path. 




1.3.3 Stationary Phase Approximation: 

The kernel of the propagator (jl.3(J|) is complex exponential of the type 
(exp(^)). The stationary phase approximation is an asymptotic ap- 
proximation of the Feynman Path Integral where h is considered as a small 
parameter. 

Consider a simple example, i.e. the evaluation of following integral 

(1.32) 

When h 0, the exponential becomes highly oscillatory function of x and 
the integral becomes nearly zero. However there are points where the oscilla- 
tions stop. These points are located where the variation of f{x) slow down. 
Therefore phase become stationary at these points, i.e. f'{xc) = 0. The set of 
points {xc} are called set of stationary points. The idea is to find stationary 
points in the domain of integration, including boundaries, and calculate their 
contribution to the integral separately. 

This is done by expanding exp {^"^-^^ around each of the stationary point 
Xc- Then equation ()1.32|1 can be written as 



1 °° 
F{—)= J dx exp 



fixc) + ^^{x-xcr+ 



:i.33) 



In most of the examples f"{x) ^ 0, and gives dominant contribution hence, 
one can neglect cubic or higher order terms. If integrand is analytically con- 
tinued into the complex plane of x, one needs to consider complex critical 
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points. Integral can then be performed by a steepest-descent method. Tun- 
nehng is one of the physical phenomenon in which complex critical points 
are important. 

By considering only quadratic terms equation ()1.33|) can be transformed 
into an imaginary Gaussian integral, which can be evaluated to give 




1 , / 2T{hi „ / if{xc) 



Turning back to equations ()1.30|) and ()1.3ip . the stationary points (or 
paths) are solutions of equation 

SR{q{t)) = 0, 

which therefore are the classical trajectories satisfying boundary condition 
q(0) = qo,q(^) = q, as well as Newton's equations. We can expand action 
W around the critical path (i.e. classical orbit) q'^ (t) 

R[cC + 5q] = R[q^] + _[q^]5q+^[q^]5q5q + 

It may be recalled that the type of extremum for the classical orbit and 
hence linear stability of the orbit depends on the nature of second variation. 
Thus the information about the stability of the classical orbit enters into the 
semi-classical framework via second variation of the action functional. The 
second variation of the action is given by 

t 

5^R[f] = J dt6ci' it) (£[/]) 5q(t); 5q(0) =6q{t) = (1.35) 



where C = —■^5q—dggV[f]6q. A solution of S'^R[f] = is a Jacobi field 
along the classical path. The equation of motion for the variation under 
this Lagrangian C is of the form of Jacobi-Hill equation. With stipulated 
boundary conditions, we have a Sturm-Liouville problem for the operator 
D= — , over the time interval (0,t). The operator D is real sym- 

metric so that it has real eigenvalues /i„ corresponding to real eigenfunction 
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u„, forming a complete basis on which it is possible to expand any variation 
satisfying boundary conditions as 



n=0 



;i.36) 



with D-u„ = /XnU„, u„(0) = u(t) = , and J^dtu^it) ■ u„(t) 
second variation becomes diagonal in this new basis. 

^ n=l 



6nm- The 



;i.37) 



The nature of the classical path depends on the sign of the quadratic form 
of the second variation, which in turn depends on the number of negative 
eigenvalues of the operator D. In general variation of eigenvalues is either 
negative or zero for solutions of equation (jl.36|) . Therefore with increasing 
time interval /i„ will decrease crossing zero at time say T„. The equation 
D-5q =0 admits a nontrivial solutions at this times satisfying b.c. (5q(T„) = 0. 
At other times no such solution exist. The conjugate points corresponds to 
the times T = T„. 

Both sides of equation p.3(i|l are then differentiated with respect to T 
and the integral ()1.35j) can be solved. 

There exist several classical paths q/ (t) that goes from qo to q during 
time t, and each of these stationary solution contribute to the propagator. 
From above comments we can write equation ()1.30|) as 

Nf 



;i.38) 



The higher order terms are obtained in ^H] • The matrix of the second deriva- 
tives ID = d'^^R[q\ is a (A^ — 1)/ x (N — 1)/ matrix given by 



D = — 



1 

At 



Ai -1 

-1 A2 -1 

-1 A3 










-1 A 







N-1 



where An are / x / matrices given by An = AajsinAt) = 26, 
the hmit At ^ the matrix D is related to the Jacobi-Hill operator D. The 



At^dlpV. In 
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propagator can be obtained evaluating all the Gaussian integrals and their 
moments. Finally one gets 



det^^ 



1/2 



exp 



^i?z(qoq;t)-^j 



(1.39) 

where ui is the Morse index. 

As mentioned earlier a Fourier transform establishes the bridge between 
the energy and time domain, the energy Green function can be obtained by 
taking Fourier transform of equation ()1.39p : 



^(q, qo) 



1 

ih 



, iEt 

dtexp(— )K(q, qo;t). 

h 



;i.4o) 



This integral can again be evaluated using stationary phase approximation. 
Here the stationary phase condition is dt[Et + Ri] = 0. If we define the 
reduced action ^^(qo, q;-E) = Et + Ri{qo, q; t(qo, q; E)) = J^^^ p dq we obtain 
a quantity independent of time, dtSi = 0. This stationary phase condition 
picks up all the classical trajectories I from qo to q at given energy E. The 
following steps are same as above and as a result we get 



1 



where 



D 



|det D\^^'^ exp 



iSi{qo,q;E) _ 

n 2 



;i.4i) 



dqodq 
dEdq 



dgodE 



Finally, to obtain semi-classical density of states as in equation (jl.23p . we 
take trace of above Green function. The contributions to the trace mainly 
comes from two sources: (i)very short paths for which the propagator is 
delta function [20] • In fact these are equilibrium points of the system, (ii) 
periodic orbits of non-zero length. The contribution of paths that are closed 
but not periodic, is negligible as a result of destructive interference among 
themselves. The former leads to Thomas-Fermi term for the average density 
of states 

da,{E) = dpdqSiE - Hip, q)). (1.42) 
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The second source leads to the oscillatory contribution for the density of 
states, which in general form can be written as 

dUE) = E E A.-(^) cos (r {Si{E) - ai)) (1.43) 

I r 

where l,r denotes the primitive periodic orbits and their repetitions respec- 
tively, and Si,ai denotes action and phase. The amplitudes Ai^r depend on 
two aspects of the periodic orbit-whether they are stable and whether they 
are isolated. The total density of states is sum of these contributions, 

d{E)=da.iE)+dosciE). (1.44) 



1.4 Density of States for Pseudo-integrable 
Billiards 

In the pseudo-integrable billiards almost all (in sense of Lebesgue measure) 
periodic orbits are marginally stable and non-isolated (i.e. occur in bands.). 
The density of states is given by [21] 

d{E) = -^ReV] [ [ dqH^^'^ (|V2mE ) exp(zai7r). (1.45) 

where m is mass of the particle (billiard ball) and I is length of the periodic 
orbits. The index aj(half of Maslov indices v) represents number of specular 
reflections at the boundary of the billiard. And h'^^ is a Bessel Function of 
third kind (Hankel function). Since periodic orbits form continuous families 
on the invariant surface and each orbit in the family has the same value of /j, 
independent of q, integration in (jl.45j) is trivial. The oscillatory contribution 
to the density of states in the asymptotic form [Ti — * 0) is then given by 

rfo.c(i?)= (^)'^'i^-^/'E^cos(^|v^+(«.-i)), (1.46) 

where represents area of the bands of the periodic orbits. The summation 
here is over all primitive periodic orbits and their repetitions. The average 
contribution from zero length periodic orbit gives the Wyle area[22j contri- 
bution to the level density 

daviE) = j^Re{ARH^'^\o)) + higher. order. terms. ,^ 

= WM0) + h.o.t. = ^ + h.o.t. > 
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where, Aji is configuration space area of the bilhard and higher order 
terms include corrections from boundary and corners of the configuration 
space. 



Integrated density of states or mode number N{E) = Number {En < E} = 
En'S>{E- En) is then given by 



Once the complete information about classical variables such as I, A is 
obtained, equations p.46|) and ()1.47|) can be used to study various spectral 
fiuctuation properties in which we are interested. 



In this chapter, we have attempted to give brief sketch of the develop- 
ments in semi-classical techniques taken place during last few decades. In 
particular the problem begun by Selberg of finding semi-classical properties 
of a quantum particle moving on a constant negative curvature Riemann 
surface, has enjoyed a fiurry of new activity in recent years. An extensive de- 
velopment of concepts by Gutzwiller and other workers have resulted in many 
applications in laser spectroscopy, Rydberg states of atoms or molecules, elec- 
tronic semiconductor devices etc. The success for the path integral approach 
is in the fact that in all these systems periodic orbits can be classified, their 
actions computed, and the path integrals can be summed. Many techniques 
has been developed to carry out this task. The Gutzwiller trace formula, 
however, gives only leading order terms in the h expansion. Higher order 
terms also has been computed recently by Gaspard|15j. 



E 
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Chapter 2 

Billiards - Classical Dynamics 



2.1 Introduction 

Billiards are dynamical systems corresponding to the motion of a point like 
particle in a bounded domain Q, which is a compact Riemannian manifold 
with piecewise smooth boundary dQ. We assume Q C Qo, Qo being a closed 
C°° Riemannian manifold. The boundary dQ consist of a finite number of 
smooth compact C°° sub- manifolds dQi, dQ2, ■ ■ -OQr of co-dimension 1. The 
points of the boundary qedQ are singular of order n if dQ is only C" at q. 
The points which arc not singular are called regular and represented by a set 
dQ = dQ \ {q of order < oo}. 

The particle reflect from the boundary according to law of elastic reflec- 
tions. The motion between reflections corresponds to the geodesic flow Gt, 
associated to a Hamiltonian Hq with Qo ^ conflguration space. The reflec- 
tion of the particle trajectories from a singular point of the boundary is not 
well defined. Being of measure zero (in Lebesgue sense) these trajectories are 
of little importance from the point of view of ergodic theory. These trajecto- 
ries may play important role in semi-classical theories, rendering diffraction 
effects. Before turning to specific types of biUiards we are interested, few 
relevant terminologies, definitions are given below. 

2.1.1 Phase space 

The phase space of the biUiards is the set of all tangent vectors of fixed 
length, M = {x = (q, v)|q e Q, v e -5'^~^(q)} where d^diuiQ and -5'^~^(q) 



33 



is the unit sphere of dimension d — 1 over q. The reflections of x e dM — 

{x = (q, v) e M|q e dQ} where dQ — dQ\{q of order 0}, are defined by 
R : (q, v) ^(q, v— 2 (nq, v) nq). Here, nq is the unit inward normal to dQ at 
q. Let TT : M ^ (^(equivalently, x = (q, v) i — > q) be the natural projection 
on Q. Since 7r~^(q) = q^S'^~^{q), we have dimM = 2d — 1. For x eM, the 
point q = 7r(x) is said to be carrier of x. M possesses a natural involution 
sending each point x =(q, v) e M into a point x' = (q, — v) e M. 

Definition 2 The measure /i on M is defined as 



where c?p(q) is the element of volume generated by the Riemannian metric 
into Q and (ia;(q) is the Lebesgue measure into S'^'^{q) 3 v. 

The flow Gt of the billiard corresponds to a vector field {G't(x), x G Mq}, 
where G't(x) is a tangent vector to Mq at x. (Mq being unit tangent bundle 
over Qq. The fiow Gt thus determines the motion of a particle with unit 
velocity along geodesic fines. We define fiow ^^at the moment of a refiection 
from boundary dQ as follows 



where, M_ = {x = (q, v) e dM \ (nq, v) < 0} and t(x) is the nearest 
strictly positive moment of a boundary refiection of the trajectory issued 
from X. 

Let Nij be the set of all interior points x gM such that the segment of 
the geodesic line in the direction of x intersect dQ on dQi fl dQj. Denote 
by N'^^\ the set of all points x eM which will be contained in VJi^jNij at 
some step of construction of the geodesic flow. Denote by N^'^\ a set all x 
for which process of construction of a geodesic leads to an infinite number of 
refiection in finite time. Then 

Definition 3 If for almost every (in the sense of measure //j x eM' — M\ 

[N^^'^ U A^*^^)), we have a geodesic segment of a finite length with end point 
located at a regular point, the billiards are said to be proper. 



djji — dp{q)duj{q) 



(2.1) 




X ^ M 

X e M 



(2.2) 
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We will consider here proper billiards only. Dynamics of billiards is of 
Hamiltonian nature. However, due to the reflections on the boundary its 
dynamics cannot be completely described by a Hamiltonian, unlike a par- 
ticle moving under the influence of a conservative force. There are many 
examples of dynamical system that can equivalently be described in a bil- 
liard. For example system of two particles of masses mi,m2 moving in a 
unit interval [0, 1], bouncing off from boundary as well as from each other 
elastically is equivalent to a triangular billiard. Another example is of three 
hard rods sliding along a frictionless ring and making elastic collisions which 
is also equivalent to a billiard on triangular table|lj. One can also deform a 
Hamiltonian system mathematically into a billiard. [21 El E] 

The billiards can exhibit all features of dynamical systems, from integra- 
bility to chaotic behaviour. The behaviour of course depends on the geometry 
of the billiard table. Billiards can be classified according to their behavioural 
pattern as: 1) Hyperbolic billiards or dispersing billiards, e.g., polygons with 
smooth obstaclesjS], some billiards with convex boundaryj^j, 2)Elliptical bil- 
liards with strictly positive curvature convex tablesJT]. 3) Parabolic billiards, 
e.g., polygonal billiards. Here we shall concentrate on polygonal billiards in 
the Euclidean plane. 

2.2 Polygonal billiards 

Let P be closed, connected, non-self intersecting polygon in the Euclidean 
plane, P C R^, whose boundary dP consist of a finite number of line seg- 
ments(edges), we denote them in arranged order dPi, dP2 ■ --dPr (some times 
we will use notations 01,02. .etc. for convenience) such that dPi±i are neigh- 
bour of dPi. The points dPi fl dPi±i are vertices of P. We thus have a set 

of regular points dP = dP = {q edP\q ^ dPi fl dPi±i}. It is obvious that 
q EdPi n dPi±i are zero order singular points. We denote its phase space by 
M (P). At each regular point q G dPi, the unit normal vector is same and 
shall be denoted by rij. 

Consider the one-parameter group {T*} of transformations on M (for 
X = (q, v) G M', — 00 < t < 00) defined as follows: T*Xo is obtained by 
starting at qo and drawing a continuous path inside P consisting of straight 
line segments and of total length \t\ and ending at x^ = (qt, Vj). The straight 
line segments (except first and last) begin and end on dP and is called as the 
link of trajectory a given trajectory. The direction change at the boundary in 
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passing from one link to the next is made in accordance with laws of elastic 
reflections. This is done via a map : — > at each point x =(q, v), q G 
dPi which acts according to the formula (TjV = v— 2 (rij, v) n. If path hits 
vertex we stop trajectory there. The set of the form 7r({T*x : — oo < t < 
oo,x G M'}) represents configuration trajectory of P. If trajectories hit 
vertices of P both in the future and in the past then trajectory is finite, we 
call it a generalized diagonal of P. 

The procedure of unfolding polygons can be applied now (see Fig. 2.1). 
The configurational trajectory vrdT^x}) has vertices on dPi. Starting from 
any reference point (qo,vo) on the trajectory, let trajectory has successive 
vertices on the boundary segments of ai, 02, 03. ..of P. Then reflections of 
P with respect to these faces transform broken line to the trajectory into 

a straight line intersecting the polygons P, Pai, -Pai,a2 Here, Pai,a2...ai is 

the polygon obtained by reflecting P with respect to the sides 
The velocity of the motion on the part of trajectory after k*^ reflection is 
Vfc = (c"afcCrafe_i •••CTaJvo- Let Gp be the sub-group of the isometry group of 
the unit circle generated by the isometries aia2..-(Jr. In the analysis of 
billiard dynamics, the singularities produced by the vertices play a major 
role. The vertex angles and the relative lengths of the edges are important 
characteristics of P as far as dynamics is concerned. If all the angles of P 
are commensurable with vr then we call P a rational billiard. The group 
Gp is finite for rational billiards. The problem of billiard dynamics deals 
with behaviour of the billiard trajectory. The question mainly falls into 
two categories, one concerning statistics which belong to ergodic theory and 
other concerning the topology of trajectories. Here we are more interested 
in the topological properties then the ergodic properties of the polygonal 
billiards. Most of the information we need for semi-classical study can be 
extracted from techniques developed in these areas. As shown in the previous 
chapter semi-classical properties of a system depends mainly on the periodic 
orbits. We shall therefore simply state some well known facts about ergodic 
properties of polygonal billiards. And discuss topological properties that 
mainly concerns about periodic orbits of polygonal billiards. 

The billiard in a typical polygon is ergodic jH]. A prevailing opinion in 
the mathematical community is that polygonal billiards are never mixing, 
but this has not been established yet. However there is a conjecture about 
existence of weakly mixing polygonal billiards jH] which is stronger property 
than ergodicity but weaker than mixing. The rational polygonal billiards are 
however, proved to be not mixing . Some of the important facts relevant to 
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a 



Fig. 2.1: Unfolding of a polygon P. Boundary segments 
are denoted by nos. 1,2, 3, 4- A trajectory starting from 
'a' in P hits sides 4-3-1-4 at b-c-d-e respectively. . 
A straightened trajectory a-b-c'-d'-e' is obtained by 
reflecting P in side 4 to get^j and so on. 



our work are presented below. 

Proposition 1 The set of generalized diagonals of P is countable.il 

Proof: Let X be the set of polygons obtained by unfolding along all 
billiard trajectories in P. Then X belongs to the set of polygons {gP : g G 
Gp} which is countable because the group Gp has finite number of generators 
and hence countable. This implies that set of pairs of vertices of any two 
polygons in X is also countable.* 

Theorem 3 For any q E P C and almost all (w.r.t. Lebesgue measure) 
\ E the closer of the configurational trajectory of the point li. = (q, v) with 
respect to the billiards flow {T*} contains at least one vertex of the polygon. 

Let 7 be periodic trajectory (of length L)m P with m links 70, 71, ■ ■ ■'jm-i 
such that 7o is the reflection of 7m- 1- Unfolding of 7 then results in set of 
polygons Pq = P, Pi, - ■ ■Pm-i.Pm- Choose a link 70 and let I be the line 
through 7o. The element g e Gp that moves P into Pm-, also moves 70 into 
7m {tilde here represents link in the unfolded polygon that corresponds to the 
link in original polygon) which belongs to the same line Z as 79. For example 
in Fig 2.1, 7 is abed.., links 70,71... are simply segments ah, he... so on and 
links in the unfolded polygons are 6c', c'd' etc. 

Thus g preserves the line /, hence g is either a parallel translation along 
/ preserving orientation, or a sliding reflection with axis / reversing the ori- 
entation. Since gP = Pm, g is the product of m reflections, m is even if g 
preserves orientation and odd if reverses orientation. 

Corollary 1 If m is even, 7 extends to a band of periodic trajectories of 
length L parallel to 7. Both boundaries are unions of the generalized diag- 
onals. If m is odd, every trajectory 7' starting close to 7 and parallel to it 
comes back after m reflections to the same edge, at the same distance from 
7 and in the same direction but on the opposite side of'j. 

Proof: It is convenient to denote 7 by 70. unfolding P along 70 we obtain: 
l)the sequence P = Po,Pi, ■ ■ -yPm, 2)the line Iq and 3) the motion g such 
that Pm = gP. Let the point xq on the edge a of P be the starting point of 
7o and let yo = gxo be the corresponding point on the edge h = ga of Pm- 
The periodicity of 70 implies that Iq goes from xq to yo. 
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Let 7 be the trajectory starting ai x E a close to Xq and parallel to 70. 
If 7 is closes enough to 70, unfolding along 7 we obtain the line / through 
X parallel to Iq and passing by the same sequence P = Pq, Pi, ■ ■ ■, Pm-i of 
polygons. Let / intersect b at y. 

If 7o is even, by previous discussion, the quadrangle formed by a, b, Iq, 
and I is a parallelogram, thus y — gx, i.e. 7 comes back to x in the same 
direction. Hence, 7 is periodic. Now start moving 70 to the right parallel to 
itself. We then get ^t, t > 0. For small t, 7t is periodic of the same length 
as 7o. Increasing t we come to the moment t = ti such that 71 = ■jn, hits 
the vertex of P. Since 7^^ is the limit of the periodic trajectory, it is also a 
generahzed diagonal or union of such. The same argument works for i < 0. 

When 7o is odd, the quadrangle above is a trapezoid and the point Xi = 
g~^y on the side a where 7 returns after m reflections is symmetric to x with 
respect to xq.» 

Theorem 4 In any rational polygonal billiard there is at most a finite num- 
ber of prime periodic orbits of odd periods. 

Consider convex rational polygon P with r vertices. Its interior an- 
gles ai,a2,- ■ -jCKr being commensurable with tt, can be written in form 
C(i = {ki/n)7r, where n > 1 and the greatest common divisor of the num- 
bers n, kik2, ■ ■ -kr is 1. Then we have following lemma for group Gp. 

Lemma 1 The group Gp is isomorphic to the symmetry group of a regular 
polygon with n vertices. 

Proof: Note that in the velocity space the composition of the sym- 
metries corresponding to the sides bounding the angle is a rotation R2ai 
of the circle by an angle 2ai = {2k.i/n)'K. Hence, for any family of 
integers sq, Si, ■ ■ ■, Sr, the group Gp contains the rotation by an angle a = 
{sqh + siki ■ ■ ■ -\-Srkr) {2TT / n) . But it follows from our assumptions that there 
exist a family sq, si, ■ ■ •, Sr, for which squ + siki • ■ ■ -\-Srkr — 1. Hence the 
group Gp contains the rotation R2-n-/n and all rotations by angles which are 
multiples of 271 /n. Moreover, the group Gp contains n symmetries with re- 
spect to the axes of the form R2TT/nl, k = 0, 1, — 1, where / is any of 
the axes of symmetry of cti, (T2, • • •, cr^. The transformations indicated above 
generate the symmetry group of a regular polygon of n sides.* 

The group Gp is thus isomorphic to the dihedral group Dn. The group 
Gp therefore has 2n elements and the circle < ^ < 27r is divided by the 
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action of Gp into 2n intervals 7r(i — l)/n < 6 < vri/n, i = 1,- ■ ■, 2n. Every 
6*, < 6* < 27r is equivalent by action of Gp to a unique 9i, < 9 < n/n, 
so that set of orbits of Gp is parametrized by [0,7r/?T,]. The orbit of every 
6, < 6 < ir/n has 2n elements, the orbits of ^ = and 9 = tt/ti have n 
elements each. 

2.3 Invariant Surface of Polygonal Billiards 

Theorem 5 If Gp is a finite group, then the billiards in the polygon P are 
not ergodic. Moreover, to each orbit of the natural action of the group Gp 
on (i.e., to the set fl = fl{yo) = {d^o ^ '■ 9 ^ Clp}), corresponds to 
the set Aq, invariant with respect to {T*}, consisting of all x = (q, v) e M 
such that \ eQ. MJI 

Proof: Suppose x = (q, v) G y4o,i.e.,v = (^qVo, Qq G Gp. Then any 
t G K, we have (qt, vj = T*x G Aq,. Since the group Gp is finite we can find 
a set of orbits of the group Gp : C C /Gp whose measure differs from zero 
or one. The set A = Uogc^c is invariant w.r.t. {T*} and yu(A) differs from 
zero or one.» 

This theorem states that for a finite group Gp, only a finite number of 
directions may be obtained when we move along billiards trajectories from 
the given initial direction. 

It may be noted that any vector v G 5^ can equivalently expressed in 
terms of angle 6. From the discussion of previous section for any 6, < 
9 < 'k/tl, the set Re = Aq is invariant under the billiard fiow. The invariant 
surface Rq are level surfaces of the function ip = M{P) — >• [0, n/n], defined as 
ijj{x,r]) = 6, where < 6' < vr/n and rj = (7^ [1011111 1121 CB]- The function is 
independent of the Hamiltonian H that gives rise to the fiow of the billiard 
and is a constant of motion. 

2.3.1 Construction of the invariant surface 

The invariant surfaces Rg for < 6 < 71 /n are isomorphic to a surface R 
which can be geometrically constructed the polygon P (see Fig. 2.2 for an 
example). To construct R which is made up of 2n copies of P, say Pi,P2 ■ 
■■, P2n, choose 9i, < 9i < n/n and denote 6*2, 6*3, ■ ■ ■92n the elements of the 
Gp-orbit of 9i in the natural order. Now take polygon Pi remove it from list 
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of 2n copies of Pi and consider any trajectory in direction 6i, it hits a side 
say a,ij of Pi. Reflect Pi on the plane about aj^, join the reflected Pi(say 
Py) to the figure under construction (i.e. R) resulting in a combined figure 
Pi + P^i = Pi and remove a copy of P^' from set of 2n copies of P. Repeat 
the whole process with P^, for trajectory in the same direction 9i but not 
necessarily along the same line. Each time trajectory meets a side Ojj, carry 
out reflection. If Pj/ is still in the list of 2n polygons join it to the figure under 
construction and remove it from the list. If Pj/ is not in the list then identify 
side Gij and corresponding side of Pj/ in the P^ which are already drawn and 
transfer a trajectory to Pj/. Repeat the procedure until all the 2n polygons 
are exhausted from the list, and resulting in a polygonal surface R. In the 
case of overlapping between P- we consider them belonging to the different 
copies of the plane. There are many possible surfaces R for a given P. The 
polygon R has an even number of sides and they are divided in pairs, each side 
of any pair differs from the other one by a parallel translation. The billiard 
flow have simple realization on R. Consider a trajectory in any direction 9, 
suppose it meets boundary segment a of R. Since sides of R comes in a pair, 
let b be pair of a. The trajectory meeting a gets instantly transferred to b 
by the parallel translation that identifies a and b. The trajectory then starts 
anew in the same direction 9 from the side b. 

If each identified pair glued together one gets the invariant surface R{ we 
keep here same notation R for simplicity reason). It is therefore clear that 
P is a closed surface without any boundary. It can also be shown that R 
is orientable. The gluing or identifying sides preserve the orientation of the 
side. 

2.3.2 Topology of the invariant Surface 

The topological type of a closed orientable surface R is determined by its 
genus g{R). The surface of genus g looks like a pretzel with g holes or 
equivalently a sphere with g handles. 

Again, let nrrii/ni, i = 1, ...,t, be the vertex angles between the sides of 
P. From discussion above one can see that each vertex aj of P with angle 
Tinii/ni gives rise to n/ui singular points aij of the flow on R. Each a^j has 2mi 
equally spaced throngs and the flow lines in a small neighbourhood of such 
singular point are shown in Fig. (2.3). Singular points of this type are called 
multisaddles, the flow lines coming into and going from the singular point 
are called the incoming and the outgoing separatrices. Thus a vertex gives 
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(A) 



(B) 



Fig. 2.2: (A) and (B) shows two different 
surfaces 'R' of a f shaped bilUard Pj . 
Here since 'n' =2, we require four copies 
ofPj, i.e. P^P^P^ and J^. Respective 
identification of the sides are shown by 
nos. 1,2... 
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rise to n/rii multisaddles on R with mj number of incoming and outgoing 
separatrices each. Varying 6 does not change the position of multisaddles Ojj 
but uniformly rotates the separatrices around aij. The index of a multisaddle 
a with 2m separatrices is equal to m — 1. The index formula for the Euler 
number gives 

' in 

x{R)=J2-a-m.)=2-2g{r) (2.3) 

1=1 * 

.and hence the genus of the surface R is 

m = i + ^t^- (2.4) 

i=l * 

Thus the topology of R is thus determined by the angles of P. The surfaces 
Rq, Rn/n are called exceptional invariant surfaces. Their topology, is not 
determined by the angles of Pjl4j. 

The topology of the invariant surface determines whether the system is 
completely integrable or not. We will now discuss this in a dynamical sense. 



2.4 Integrability and beyond: 

The dynamics of classical Hamiltonian system is completely described by the 
Hamilton's equations of motion 

dH , dH 

Qi = — ; and pi = - — 
opi oqi 

where symbols have their usual meaning i.e. H is the Hamiltonian and (gj, Pi) 
are the canonical coordinate-momentum pair. 

Consider a system of / degrees of freedom with the Hamiltonian H. The 
Hamiltonian induces the flow on the phase space M. The dynamical sys- 
tem is integrable (in the sense of Liouville) fSl ^] if there exists a canonical 
transformation such that the new momenta are the constants of motion. A 
necessary and sufficient condition for this transformation to exist is that it 
should be possible to find / analytic functions {Fi}^^^ ^ on M satisfying 
following conditions : 

• The Poisson bracket, {Fi,H} = for all i, which implies Fi{^tx) = 
Fix) for all a; e M, 
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Fig. 2. 3 Vector field near singular point with m 
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all F are in involution with each other i.e. {Fj, Fj} = for all i,j = 
1-/, 



• all grad{Fi) are linearly independent. 

/ 

Here the Poisson bracket {A, B} = - ^ff), and grad = 

(Vg,Vp). One can then define a level set of the functions {Fi} as = 
{x G M\Fi{x) = ki^ i = 1.../}. This level set is invariant under fiow $f. fur- 
thermore if Mfc is compact and connected ,then it is diffiomorphic to the 
/—dimensional torus [IHj. The vector field {Vj} ; Vj = (VpFj, — V^Fj) is 
linearly independent, tangent to M^, and commutating. According to the 
Noether's theorem [T71 HH] these constants of motion (or isolating integrals) 
result from the symmetries of the dynamical system. 

Polygonal billiards are examples of dynamical systems where there exist 
one or two constants of motion depending upon the angles and ratio of the 
sides of the polygon. For some rational polygons with a rational ratio of 
their sides, there exist two constants of motion satisfying all integrability 
conditions except one i.e. {Fi,Fj} 7^ at the vertices of the polygon. Thus 
the vector field V becomes singular at countable number of points in the 
phase space. Hence it is not possible to obtain global action-angle variables 
or constants of motions. Following Richens and BerrvfT^. we call these 
systems a pseudo-integrable. 

2.5 Birkhof-Poincare Maps and Interval Ex- 
change Maps 

The Birkhof-Poincare map is first return map for billiards. For simplicity 
let P be simply connected billiard. We enumerate the vertices of P coun- 
terclockwise Ao, ...An-i, An = Aq. Denote by aj = [y4j_i, Aj], i = 
the edges of P. The boundary dP with length coordinate x is isomor- 
phic to the circle of perimeter equal to the length(L) of dP. The set Q C 
M{P) of vectors with footpoints in dP can be parameterized by coordi- 
nates < X < L and — 7r/2 < 9 < 7i/2. The vector {x,6) has foot points 
in a; G flj C dP and angle 6 between it and the inner normal to a. Set 
Lk = + ■ ■ ■ + \ak\ , k = 0, ...,n. The coordinates 9 in {x,9) is not well 
defined for x = Lk because these x are vertices of P and the angle 9 can be 
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measured with respect to any of the two normals. The set with points ex- 
cluded is isomorphic to the cylinder [0, L) x [— vr/2, 7i/2] with deleted intervals 

Lox[-7r/2,7r/2], , L„_i x [-7r/2, 7r/2] (see Fig. 2.4). The Birkhof-Poincare 

map F : Q Q is defined as follows: consider a trajectory from x in direc- 
tion 6, [x, 6) G fl. When trajectory hits dP the first time and bounces off, it 
determines another point {y, r]) = F{x, 6) to Q. The mapping F is not well 

defined on the deleted intervals Lq x [— 7r/2, 7r/2], ,Ln^i x [—h/2,tt/2] as 

well as on the boundaries 6 = ±7r/2. Fix an edge = [^i-i, Ai] and a vertex 
Aj^ j ^ i — Points (x, 9), x E Qi, such that the ball goes to the corner 
Aj form a curve in Q on which map F is not well defined. Each rectangle of 
Q is divided by these curves into the domains of continuity of F. Thus, the 
set of discontinuities of F is the union of a finite number of curves in Q, F 
is obviously invertible and F~^ is the Birkhof-Poincare map for the billiards 
with time reversed. The F -invariant Lebesgue measure is cosOdOdx. The 
properties of the billiard fiow are easily translated into the properties of the 
mapping F. 

2.5.1 Interval exchange Map For Rational Billiards 

Since the billiard fiow G* decomposes into the family Gg of fiows, < < 
vr/A^, the Birkhof-Poincare map F also decompose into the one parameter 
family Fq : VLq ^ VLq mappings, where fig = i7 fl -Rg is set of vectors with foot 
points on dP with directions equivalent to 9 and Fq = F \ VLq. The family 
{Fg} is the family of interval exchanges 12011^ ■ 

Suppose the space / is semi-interval [0, 1) and ^ = {Ii....Ir) is a partition 
of / into r > 2 disjoint semi-intervals Ji = [0, di), = [dr-i, 1), < di < 
... < dr^i < 1. Let u = (tui, ...Ur) be permutation of the number {1, ...r}. 

Definition 4 Suppose the transformation T : I —>■ I is a translation T^.x = 
X + aj(mod l)on each of the semi-intervals Ii(the number ai depends on i) 
and exchanges the semi-intervals according to permutations u i.e. The semi- 
intervals T li = Taji = /j adhere to each other in the order : then T 
is said to be interval exchange transformation corresponding to the partition 
^ and the permutation ut. 

Thus if / is piece of wire then the transformation T essentially is the 
cutting of I into pieces Ii....Ir, rearranging them according to permutation lj 
and welding them together again. Identifying / with the circle we can think 
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6l=7r/2 



6l=-7r/2 





1 1 1 
1 1 1 

! I 1 

1 1 
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1 ! 1 

—•--4--^ 

1 ! ! 




1 

92 \ 1 1 




1 I 
IJ., 1 

1 1 1 
1 1 1 





Fig. 2.4: Domain ofBirkhof-Poincare mapF . 
Vertical dashed lines (corresponding to x=Lo... x 
are deleted intervals from the domain. See text 
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of T as an interval exchange on the circle. An exchange of two intervals is 
given be one number < a < 1 and it is simply rotation by angle a. It is 
clear that interval exchanges are invertible transformations of / preserving 
the Lebesgue measure p and the numbers a\....ar are well defined (mod 1 ) 
by the pair oj). 

Now to see that the family \Fq\ is the family of interval exchanges recall 
that in coordinates (x, ^) the set f2 is [0, V) x (— 7r/2, 7r/2) and F{x, 6) = {y, rj) 
where t] locally depends only on 6. The map F perseveres dm = sin 6dx which 
is interpreted as the mass element carried by the flow in direction 6. Fix a 
direction 9 and for each side of P let —tt/2 < < ... < Oi^^ < n/2 he the 
directions DAr-equivalent to 9 (they depend on i because Oj determines angle 
of reference). The set fig is the union ai x {9i^...6ij^} U ... U Op x {9p-^...9pj^} 
of horizontal intervals. The mapping F preserves the set Qg which is union 
of pN horizontal intervals (they can be glued into one interval) and the 
length element dm on Qg. The positive orientation of P induces orientation 
of Qg (each interval oriented from left to right), the map F reverses the 
orientation. Thus, the restriction Fg of F to Qg is an interval exchange with 
the flipping of intervals. Multiplying Fg by the trivial orientation reversing 
map J : (x, 9) ^ {L — x, 9) we obtain honest interval exchange. Dividing by 
total length mg of Qg which is 

mg = \ai \ (sin^i^ + ... + sin6'ij^) + ... + \ap\ (sin^p^ + ... + sm9pj^) 

we normalize Fg to an interval exchange on [0, 1) (with flipping). For a fixed 
polygon P the parameter of Fg i.e. the number of exchanged intervals, the 
permutation u and the length of interval depends only on 9. The obvious 
upperbound on the number of exchanged intervals is p'^N. 

The interval exchange maps has been used to prove many statistical and 
topological properties of polygonal billiards|2ni 1^- Our main interest is in 
enumeration, classification of the periodic orbits. In the next chapter we will 
use modified form of interval exchange to carry out this task in case of some 
specific examples. 
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Chapter 3 



Periodic Orbits in Some 
Pseudo-integrable Billiards 

3.1 Introduction 

In this chapter we take up first task to study spectral fluctuations of pseudo- 
integrable billiards within semi-classical framework, i.e. to enumerate and 
classify all the periodic orbits of some typical pseudo-integrable billiards and 
also define their respective actions. We shall consider 7r/3-rhombus biUiard 
for this purpose. This specific example serves as a paradigm model as there 
are hardly any general results known in the literature. Furthermore the 
methodology we develop here, being based on interval exchange transforma- 
tion, may be used to enumerate and classify periodic orbits in many rational 
polygonal billiards. 

3.2 Periodic orbits of the 7r/3-rhombus bil- 
liard 

The 7r/3-rhombus biUiard is an example of an almost integrable system. As 
discussed before, billiard flow for a particle inside an almost integrable poly- 
gon is called an almost integrable billiard. The most fascinating mathemati- 
cal questions are related to the periodic orbits of these billiards. For complete 
analytical semi-classical study of such systems one needs to enumerate and 
classify all the periodic orbits and also able to compute actions of periodic 
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orbits of the given system. The example where the enumeration and classifi- 
cation of periodic orbits is analytically carried out are very rare. Apart from 
the trivial enumeration of the orbits for a separable barrier billiard there is 
no instance where a complete study exists. 

3.2.1 Enumeration 

To begin with, let us briefly recapture how the motion of a particle inside 
7r/3-rhombus shaped encloser can be visualized as motion on an equivalent 
barrier billiard It is simple to see (ref. Fig. 3.1) that after three suc- 

cessive reflections {A B ^ C) of the rhombus A around a vertex of angle 
27r/3, the rhombus returns upon itself but with reversed orientation with ver- 
tices 2 and 4 interchanged. In other words, we obtain the final configuration 
of vertices as if we have reflected the rhombus about the shorter diagonal 
of the rhombus joining the vertices 1 and 3. If we continue the reflections, 
it will take exactly three more, or equivalently, another reflection about the 
shorter diagonal, for the rhombus to identify itself with original orientation. 
In this picture, due to double-valuedness of configuration of vertices per di- 
rection ( by direction we mean one of the three directions the rhombus is 
facing in Fig. 3.1), one can visualize three rhombus-orientations in Fig. 3.1 
on one sheet (or, plane) and the subsequent three orientations (required to ob- 
tain the original configuration of vertices in rhombus-A) on another sheet (or, 
plane). One can visualize a trajectory of particle reflecting from a wall of the 
rhombus by letting the particle move straight and appropriately reflecting 
the rhombus about the wall. It is this way of analysing that turns out to be 
more fruitful and hence the discussion on the tessellation of plane by rhombi. 
Due to its equivalence to the Riemann surface of z^/^, z being a complex vari- 
able, we notice that the two sheets discussed above are joined along straight 
lines (The complex counterparts are branch cuts) that cannot be crossed; we 
call these barriers. Furthermore, as we have seen above, going to the next 
plane is to compensate for a phase tt, the trajectory must reflect from the 
barrier. Alternatively, after three reflections, we can reflect third rhombus 
back onto itself (hence compensating for the phase in the third step), i.e. 
the fourth rhombus comes to lie under the third rhombus. If we continue 
reflections now, the sixth rhombus will come to lie under the first rhombus. 
In this picture the point in Fig. 3.2 will become a (monkey) saddle point. 
Continuing the process of reflections, we can get planes connected by cuts 
between the saddle points. This construction, projected onto two dimen- 
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Fig. 3.1: Refection of the rliombus about a 2Ji/3-vertex. 
The vertices 2 and 4 are interchanged on third 
reflection. 
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Fig3.2: Apearance of a saddlepoint a the vertex 1, 
that results in the cuts in fig. 3.3 
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sions, entails an orbit looking like a zig-zag line. We can now construct the 
fundamental region using six replicas of rhombus, and subsequently tessel- 
late the two-dimensional plane by staking the fundamental regions side by 
side, exploiting the translational symmetry. On doing so, we will generate a 
barrier billiard shown in Fig 3.3 with two sets of planes (call them top and 
bottom) interspersed with each other. 

Now concentrate our attention on the barrier billiard where the barrier 
to gap ratio is two. the barrier constituted of contributions from two rhombi 
and hence, there are two distinguishable sub-barriers giving rise to a single 
barrier of length twice that of gap. Classification of the periodic orbits had 
been carried out for a barrier billiard with barrier to gap ratio equal to 
unity It must be noted that the barrier billiard corresponding to the 
rhombus problem is more general than the barrier billiard studied in 
since the barrier to gap ratio is two and moreover the barriers appear in an 
oblique manner at an angle of 7r/3. 

First, we observe that bifurcations of the orbits take place at the two 
ends of the barriers and at the center of the barrier (e.g. one example of such 
bifurcation is shown in Fig. 3.3, where two neighbouring trajectories on the 
opposite side of a trajectory OAIJC. bifurcate at end of a barrier near vertex 
C). This is due to the fact that each half of the barrier is contributed from 
two different rhombi in the fundamental region, and the point of bifurcation 
actually corresponds to a vertex. 

The single connected surface is made up of two planes - a top and bottom. 
Under the covering of the surface by fundamental regions (double hexagons), 
the surface divides into alternate arrays of both planes containing barriers. 
Obviously, it does not matter which plane is called top(or bottom). This 
argument allows us to choose an origin which, for obvious reasons, dictated 
by symmetry of the barriers, is chosen to be the center of the barrier, denoted 
by O in Fig. 3. 3. Calling the length of a side of the rhombus by L, the barrier 
length is 2L and gap length is L. On the vertical axis, the perpendicular 
distance between adjacent arrays of a (top/bottom) plane is y/SL. Since the 
factor of is common in the vertical axis, we choose to measure the length 
in this direction in terms of \/3, thereby making the ordered pairs labeling 
the points purely consisting of integers, {q,p)- For instance, a point O' in the 
Fig. 3.3 will be labeled by (3, 2). A typical trajectory on this surface will be 
made up of alternate motions in the top and the bottom planes via barriers 
(i.e. when a trajectory hits a barrier on a upper plane it simply appears on 
the lower plane and vice versa). Since every plane consist of an identical 
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Fig. 3.3: Tesselation of the plane by stacking the fundamental regions 
exploiting the translational symmetry. Barriers are shown by thick 
lines. Trajectory OABC... can be straighten as OA'B'C. . Two nearby 
trajectories parrallel to OABC are also shown. Owing to the symmetry 
of the barriers, the integer labelling can be done, e.g. the point 0' can 
be labelled as (3,2). 
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array of barriers, the trajectory starting at an angle with the plane from an 
initial point and ending on an equivalent point on the same plane constitutes 
a periodic orbit (see e.g. Fig. 3.4). Instead of following the zig-zag path, we 
can unfold the trajectory into an exactly equivalent straightened version as 
shown in Fig. 3.3. As stated above this straightened trajectory will lie on 
both planes, crossing planes at the barriers. Subsequently, we must decide 
which directions lead to periodic orbits. By virtue of the integer labeling, it 
is clear that all those directions that end on (integer, integer) ordered pair 
correspond to periodic orbits. Leaving apart a factor of \/3 in the vertical 
direction, these directions correspond to rational gradients on each (top, 
bottom) plane. Also, starting from origin, we must consider only those end 
points such that q and p are co-prime, since only such pairs results in primary 
periodic orbits. If q and p are not co-prime but have g.c.d. /c, then orbit 
ending at (g, p) represents A;*'' repetition of a primary periodic orbit ending at 
{q/k,p/k). By the symmetry of the barriers on the plane, we need to restrict 
ourselves to an upper-half region only. Further taking care of geometry of 
our system, we restrict further to either p < q or q < 3p, obviously then, 
three classes emerge, viz. {q,p) can be (i) (odd, odd), or, (m) (odd, even), or, 
(Mz)(even,odd). Next, we have to classify the number of bands or families of 
periodic orbits that correspond to each direction. 

For equal barrier to gap ratio, it was shown that any trajectory with irra- 
tional gradient can be approximated arbitrarily well by trajectories with ra- 
tional gradients, utilizing Klein's string construction or the continued fraction 
expansion. However, those trajectories never close in position and momen- 
tum both, rather they form a curious zig-zag path, for quadratic irrational 
gradients, has a fractal dimension. For the barrier in our case, same holds. 
Hence, we conclude that the trajectories with any irrational gradient do not 
close. It, therefore follows that if we take into account all the rational gradi- 
ents avoid over-counting and classify different families/bands, we would have 
enumerated all the periodic orbits. 

We shall now use modified form of interval exchange transformations and 
name it polar construction for the reasons of clarity and easy generalization. 

3.2.2 Polar Construction 

Exploiting the periodicity of barrier-gap-barrier-gap... string, we wrap the 
basic string of skewed sets of barrier-gap-barrier on two circles representing 
two different planes. Each circle has three basic divisions, coming out of 
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Fig. 3.4: The periodic orbit (ADCBA) in thie rtiombus 
abed is equivaient to ttie periodic orbit (ADC'B'A') 
in ttie barrier biiiard, wtiicti doses on equivaient 
point on ttie ottier plane. 



58 



two barriers(from two rhombi, joined together) and a gap. Each division 
now corresponds to an angle 27r/3. The fact that barriers on two planes 
are stacked in a skewed manner is accounted for giving appropriate phase 
difference between equivalent points on two circles (see. Fig. 3.5). A {q,p) 
direction can be represented these circles by the following procedure. 

Divide each segment of the inner (outer) circle into p parts, after fixing 
the origin at the point joining the two sub-barrier segments. The origin of the 
outer (inner) will be at an angle of 7iq/3p from the origin of the inner (outer) 
circle moving along the circular arc in a definite sense (we use anticlockwise 
movement). After fixing the origin, the outer circle has to be divided into 
3p parts of equal length with one division at the origin. Since there are two 
gradients, positive and negative, we adopt the convention of marking outside 
(inside) of the circle as a representative of positive (negative) gradient. 

To follow a trajectory, we start with an arbitrary point on one of the p 
sub-segments of the barrier segment on the right of the origin, on the in- 
ner circle. Next point will be on the outer circle just at the same distance 
(number of sub-segments) from the origin (of the outer circle) as the previous 
point was from its respective origin. As is clear from earlier discussion, these 
points must alter between the outer and the inner circles. The following point 
comes on the inner circle q segments away from the original one, and the next 
on the outer circle q segments away from the earlier one, and so on. Going 
on in this way, after a finite number of points, we will reach the starting 
point on the inner circle and that would make one periodic point. Although 
this is generally valid, in some cases it may not lead to the minimal length 
of the orbit. Such would be the case when, exactly after half the number of 
reflections, the trajectory will close, i.e., reach a point corresponding to the 
same respective sub-segment as it started with, on the other circle. At first 
sight, it might seem erroneous to consider this as a periodic orbit. However, 
it may be recalled that though the procedure of erecting barriers and subse- 
quent polar construction is a way leading to an easy classification of orbits; 
however, more fundamental idea is the straightening of a trajectory-reflecting 
the domain about the edge on which the particle is incident. For the cases 
where the length of an orbit turns out to be double by polar construction, one 
can easily see that using the domain-reflection method, one gets the correct 
length. Thus, without ambiguity, the trajectory in the polar construction 
must be considered periodic even if it seems to be closing on the other plane 
(see, for example illustration in Fig. 3.4). With this clarification, we need 
only think in terms of the polar construction. Different sequences of barriers 
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and gap correspond to the different orbits. 
3.2.3 Classification 

As we have seen above, the most elementary classification is in {q,p) being 
{odd, odd), {odd, even) and {even, odd). Having set the origin at the center 
of the barrier, the trajectory sets off in some rational direction and reaches 
either a center of some barrier or end (left or right) of some barrier. It is 
rather obvious to see that for each of the three cases written above there are 
subclasses which we shall call: the center-to-center(CC) case and the center- 
to edgc(CE) case. Our Procedure of classification is in following steps: 

(i) using the polar construction, we depict the trajectory on the circles, 
with an opening and an ending point on one of the sub-segments of the 
segment. We go on to the other to other sub- segments of the same segment, 
exploring the positive and negative gradients till both sides of the circles are 
filled. In all, we must fill 12p points. 

(ii) check if the orbit has already closed a half-way on the outer circle or, 
equivalently, on the other plane. 

With these steps in mind, we now take up each class separately and 
classify the bands of periodic orbits in full. 
Case 1: Odd-Odd CE 

We first describe through a simple example as to how we would arrive at 
general conclusions. Our approach would be to make a conclusion based on 
empirical data obtained by "brute force". At the end, we will provide with 
a rationale supporting and explaining the conclusion obtained. 

Let us consider the case of {q,p) = (1,3). The corresponding polar con- 
struction is shown in Fig.3.5. As can be seen, point 1 and point 10 identify 
with each other, forming a periodic orbit after six bounces. Also, there is an 
orbit with negative gradient. It should be noted that all the sub-segments 
are visited by just these two orbits. The orbit close a half-way of {3q, 3p),on 
the other plane and, there are two bands of orbits. Drawing the polar con- 
struction for other odd-odd CE cases, it can be seen that there are only two 
bands of orbits as seen above. 

Now we come to discuss the rationale behind this classification of orbits. 
Due to the polar construction, each sub-segment is an arc of angle 2n/3p. 
Translating the formation of periodic orbit by the polar construction into an 
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Fig. 3.5: Polar map for the odd-odd center-edge case 
(q,p) = (1,3) 
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equation, we trivially get 



A^(27r/3p)g = 27rM 



(3.1) 



where denotes the number of sub-segments and M denotes the number of 
rotations by 2tt. Henceforth, we call by the "crossing index" and M by 
the "rotation number" eq. ()3.1|) is simply 



where M,q,p are positive integers. Since this is a CE case, q is not a 
multiple of three, i.e. q ^ 3l,leZ. Thus, the only way in which eq. (j3.2|) can 
be satisfied is if g = M and 3p = N. Note that M and N will be odd as both 
q and p are odd. The crossing-index on one circle is 3p implying that the 
total crossing- index is 6p after which the orbit closed. In all, there are 12p 
sub-segments and hence there are exactly two bands of periodic orbits. 

In general, the crossing-index is given by 



where [a, b] denotes the lowest common multiple of a and b. Trivially, for 
the CE case, N = 3p and M = g,for the CC case, N = p and M = q/3. 
Case 2: Odd-Odd CC 

Firstly, all the rational directions pointing toward an infinite number of 
avenues correspond to the bands of periodic orbits. Points (g, p) correspond- 
ing to avenues are of the form (3(2fc — 1), 1), keZ^. For each direction there 
will be two bands corresponding to positive and negative gradients. 

We go over to a representative of a general case, viz., {q,p) = (3,5). 
The polar construction is depicted in Fig. 3.6. There are three strings 
of points corresponding to three bands of periodic orbits. Including the 
opposite gradients, there are four distinct bands in all. The string starting 
with point 1 and ending with point 11 closes after six reflections. Same is 
the case for IB — IIB orbit. The opposite gradient counterpart of 1 — 11 
{IB — IIB) is equivalent to itself, starting from the outer circle (at 6(65)). 
So these strings give us two bands of periodic orbits. Care must be taken 
for the periodic orbit starting with the point lA. The orbit closes at 6A 



Nq = 3Mp 



(3.2) 




(3.3) 



M = [q,3p]/3p 



(3.4) 
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as this is the subsequent corresponding to lA on the inner circle and the 
gradient matches. The orbit 1 A — 6 A closes after four reflections and occurs 
in a band. Taking the opposite gradient, we get two bands of periodic orbits 
here. For other odd-odd CC cases, we obtain the same results as in the above 
mentioned case. Let us now see the reason for the occurrence of four bands 
of periodic orbits in this case. 

From eqs (j3.3|) and (j3.4j) . M = q/3 and N = p. For each orbit, the 
crossing-index will be 2p accounting for the other circle also. In total, there 
are 12p segments and it clearly follows that there must be six bands of pe- 
riodic orbits. Subtracting the two equivalent bands, we are left with four 
bands. 

Case 3: Odd-Even CE 

Consider {q,p) = (1,2). The polar construction is drawn in Fig. 3.7. 
There is a bifurcation in the band of trajectories starting with point 1, and 
further continue the two bands, primed and unprimed, to eventually close 
at 13' and 13 respectively. This feature, which can be succinctly put as 
bifurcation of the vector fields at vertices and continuation of trajectories 
in the form of bands (a signature of zero Liapunov exponent), is typical of 
pseudo-integrable system. The points 1 and 6' are, indeed, identical. Orbits 
emanating from 1 and 6/ {lA and 6 A') will be the same . Consequently, 
on allowing the opposite gradients, we get just two bands of periodic orbits. 
All other examples of this class give rise to same number of bands and the 
orbit-types are also similar. Of course, the lengths and other details will be 
different. 

As the arguments for justification follow on the same lines, we do not 
repeat them for this and further cases. 
Case 4: Odd-Even CC 

For this case, consider {q,p) = (3,2), polar construction for the same is 
drawn in Fig. 3.8. Due to bifurcations, we have drawn double the marks, 
thus, making explicit that we have to fill 24p points in all. As can be seen from 
the diagram, there are four bands of periodic orbits: an orbit corresponding 
to the string 12. ..9; another orbit corresponding to the string 1A2A...9A] 
and a positive-negative gradient pair, 1'2'...5' (plus the opposite gradient). 
It is interesting to see that bifurcations of these vector fields take place at 
all possible places, viz. Zi, Z2and Z3. Considering other odd-even cases, one 
can easily see that the same conclusions about the number of bands etc. will 
hold, leaving apart the finer details. 

Case 5: Even-Odd CE 
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Fig. 3.6 : Polar map for the odd-odd center-center case 

(q,p) = (3,5) 
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Fig. 3.7 : Polar map for the odd-even center-edge 
case (q,p) = (1,2) 
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Fig. 3.8 : Polar map for odd-even center-center case 
(q,p) = (3,2) 
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Consider {q,p) = (2,1), the polar construction is shown in Fig. 3.9. 
There are two bands of periodic orbits. In fact, the periodic orbit formed with 
opening point as the negative-gradient-equivalent of the point 1 is identical to 
the orbit 1A2A...7A, starting from 2A. Thence, strings 1 2. ..7 and its negative 
gradient counterpart are the two bands in this case. Case 6: Even-Odd 
CC 

A general case can be studied through the example of {q,p) = (6,5), 
polar construction being depicted in Fig. 3.10. Although the diagram is 
getting rather complicated, in the same manner as discussed earlier, it can 
be concluded that there are four bands of periodic orbits. General validity of 
this conclusion can be verified without undue hardship. The results of this 
section, along with the lengths of the periodic orbits and the phase space 
areas of the bands in which they occur are summarized in Table I. This 
concludes our discussion on the enumeration and classification of the periodic 
orbits in 7r/3— rhombus billiard. The periodic orbits of other billiards such as 
Hannay-McCraw billiard |3j or rectangle billiard with slit at the center can 
also be enumerated and classified in the same way. 

3.3 The H-M Billiard 

The particle moves here in a configuration space where a hard line-segment 
barrier is placed at the center parallel to one of the sides of the periodic cell of 
side length L, and the length of the barrier is half of the side length(i.e. L/2). 
The invariant integral surface of this system has genus equal to two {g = 2). 
The fundamental region consist of two cells, due to presence of a barrier (see. 
Fig. 3.11). The opposite boundaries of the fundamental region are identified 
topologically resulting in one handle of a sphere. The barriers in each cell 
of the fundamental region has to be identified topologically which results in 
another handle of the sphere. By stacking the domains of the billiard side 
by side in both the orthogonal directions, one obtains an infinite lattice of 
barriers and gaps, with barrier to gap to ratio unity. One can label the end 
points of barriers by integer-pairs which form lattice points. It can be easily 
seen that the straightened version of a rational gradient (= \p/q\) trajectory 
will initially meet lattice point (g, p) and then repeat itself by meeting lattice 
points {mq, mp) where m ^ Z. On the other hand, the irrational gradient 
trajectory will never visit any lattice point though it will come arbitrarily 
close to many lattice points, hence will never be periodic. Thus the periodic 
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lOB' 16 



Fig. 3.10:Polar Map for the even-odd center-center case 
(q,p) = (6,5) 



69 



orbits in the system are the ones which hit any lattice point {q,p) in this 
array of barriers, the gradient of such trajectories will be given by \p/q\- By 
above arguments, we need to consider only the pairs(g,p) such that q and 
p are co-prime since they only give a primitive periodic orbits, and points 
{mq,mp), where m ^ Z, gives m repetitions of a primitive periodic orbit 
corresponding to {q,p)- Each such {q,p) gives different number of bands or 
families of periodic orbits, depending on whether the pair is odd-odd (o, o) 
or even-odd (e, o) or odd-even (o, e). The length of periodic orbit in a given 
family corresponding to a lattice point (g, p) is given by 

l = cL^{q^+p^) 

where c depends on the number of families of periodic orbits. The periodic 
orbits can now be enumerated and classified exactly same way as done in the 
case of the 7r/3-rhombus billiard |3]. The result is summarised in the table 
3.2. 

3.4 The Single Slit Rectangle billiard 

Another simple example of pseudo-integrable billiard is a simple modifica- 
tion of H-M billiard, where instead of periodic cells one can consider a linear 
barrier of length L placed at the center of the rectangle (L,2L), parallel to 
longer side. This simple modification results in the pseudo-integrable billiard 
whose invariant surface is topologically equivalent to a sphere with three han- 
dles(genus, g = 3). The fundamental region now consist of four replicas of 
the configuration space instead of two as in the H-M billiard. Opposite sides 
of the fundamental region are identified as well as each pair of barrier lying 
one above each other has to be identified separately which results in genus 3 
surface (see Fig. 3.12). By stacking the domains of the billiard side by side 
in both the orthogonal directions, one obtains an infinite lattice of barriers 
and gaps, with barrier to gap to ratio unity, which is exactly same as the 
lattice of barriers and gaps one obtains in the case of H-M billiard. One 
can proceed in the same way as before to label lattice point. The periodic 
orbits can then be enumerated and classified using polar maps. The result is 
summarised in table 3.3. In this chapter we have developed a methodology 
to enumerate and classify periodic orbits of some pseudo-integrable billiards. 
The same methodology can be used to study large number of rational polyg- 
onal billiards in particular one can apply technique to rectangular billiard 
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Fig. 3.11: H-M billiard (pqrs). Fundamental region requires two 
replicas as shown. V,H represents identified sides of a period- 
ic cells. B,B' are barriers in respective cells. A trajectory 
abcdefg... starts at 'a', meet H at 'b', due to identification it is 
then transferred to 'c' and so on. Note that 'd' is equivalent to 'e' 
since lower side of B is identified with upper side of B' and 
vice versa. 
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Fig. 3.12: The fundamental region of a single slit rectangle 
billiard which requirs four copies of configuration space A^. 
Side identification is indicated by 1,2,3,4 and barriers by B 
and B'. 
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with many slits(or barriers) or L-shaped billiards. 
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Type 


Classes 


JNo. oi 
Families 


Closing 
Point 


Length 


Band 
Areas 


Center-center 














odd-odd 


"1 1 

1-1 


(g/2,p/2) 










2 




Lfpq 






odd-even 


1-1 




Lfpq 








2 


(2g,2p) 


'^Lfpq 


'2Ar 




even-odd 


1-1 




Lfpq 


Ar 






2 


(2g,2p) 


'^Lfpq 


2Ar 


Center-edge 














0(1(1-0(1(1 


1-1 


(3(7/2. 3iV2) 


■iLf,j2 


3-4,, 




odd-even 


1-1 


(3g, 3p) 


"^Lfpq 


^Ar 




even-odd 


1-1 


(3g, 3p) 


^Lfpq 


^Ar 



Table 3.1: Summary of the results obtained in this chapter about periodic 
orbits of 7r/3-rhombus billiard. The symbol fpq = {q^ + 3p^)^/^, L and Ar — 
31/2^2^2 are length and area of the biUiard. Note that total band area 
occupied by periodic orbits assigned to a single lattice point is always 6Ar 
(i.e. the phase space area of the biUiard). 
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Classes 


No. of 
Families 


Closing 
Point 


Length 


Band 
Areas 


odd-odd 


1 


{Aq, 4p) 




2Afl 


even-odd^ 


1-1 


{q,p) 


Lfpq 


Ar/2 


even-odd^ 


1 


(2?,2p) 




Ar 


odd-even 


2 


(2g,2p) 




Ar 



Table 3.2: Summary of results for single H-M billiard, The symbol fpq — 
[q^ + p^)^/^, L and Ar — are length and area of the billiard. Again here 
total band area occupied by periodic orbits assigned to a single lattice point 
is always 2Ar (i.e. the phase space area of the billiard). 



Cleisses 


No. of 
Families 


Closing 
Point 


Length 


Band 
Areas 


odd-odd 


2 


(4g,4p) 


4L/p,/2 


2Ar 


even-odd^ 


2 


(2?,2p) 




Ar 


even-odd^ 


2 


(2g,2p) 


'^Lfpq 


Ar 


odd-even 


2 


(4g, Ap) 




IAr 



Table 3.3: Summary of results for single slit rectangle billiard, The symbol 
fvi ~ il"^ + P^y^"^, L and Ar = 2L^ are length and area of the billiard. Here 
also total band area occupied by periodic orbits assigned to a single lattice 
point is always 6^4 (i.e. the phase space area of the billiard). 
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Chapter 4 

Growth Rate of Periodic Orbits 



4.1 Introduction 

An important characteristic of a dynamical system is the asymptotics of the 
number of periodic orbits(or growth rate of periodic orbits) having lengths 
less than or equal to /. We denote this function by F{1). From the point of 
view of application, semi-classical theory of spectral statistics can be carried 
out only after knowing the exact form of this law. The problem addressing 
the distribution of the periodic orbits with the period of the orbits is discussed 
in this chapter. The answer to the problem is known only for the chaotic 
dynamical systems which possess only unstable, isolated periodic orbits. In 
such cases growth rate of periodic orbits is known to be exponential^ E]. 
e.g. For so called axiom A systems F{1) is given by 

F{1) = exp {kI} / K,l as I ^ oo 

where k > is the topological entropy of the system. For the pseudo- 
integrable and their special cases, the almost integrable billiards(A polygon 
P is called almost integrable if the group Gp generated by the reflections in 
the side of P is a discrete subgroup of the group of motions of the Euclidean 
plane R"^), there are some conjectures and incomplete results. 

For example a theorem due to [SI E] states that in case of any arbitrary 
polygon P, the function F{1) grows slower than any exponential, as / oo. 
By Corollary (2.1), F{1) is bounded above by the number of generalized diag- 
onals of length less than /. The sub-exponential growth rate of the number of 
generalised diagonals is estimated from above by the entropy of the billiard. 
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which is zero jSJ E]. The above theorem is the only known upper bound on 
the periodic orbits in general polygons. As for a lower bound, no result has 
been established yet. On the other hand, periodic orbits in rational poly- 
gons have efficient bounds from below and from above. The proofs of these 
estimates use the theory of holomorphic quadratic differentials on Riemann 
surface. We state here a theorem [7| |H] regarding this. 

Theorem 6 Let P be a rational polygon. Then there are positive constants 
< ci < C2 such that for all sufficiently large values of I, we have 

cif < F(l) < 02^. (4.1) 

These quadratic bounds are not likely to hold for general (irrational) 
polygons. The expectation is that, for general polygons there are polynomial 
bounds on F{1). 

Conjecture 1 Let P he an arbitrary polygon, then there exist positive con- 
stants Ci,C2 and integers 1 < ni < n2 such that for sufficiently large I we 
have 

an < F{1) < 02^. (4.2) 

For arbitrary rational polygons condition ()4.1|) may well give best possible 
estimates on the F{1). However, there are non-trivial examples when the 
asymptotics of F{1) can be computed exactly. 

Our investigations on some models are carried out analytically with com- 
parisons shown with the numerical results. We sharpen the existing theorem 
on proliferation law for almost-integrable systems in a significant manner. 
Our analysis enables us to give a general law of proliferation of the periodic 
orbits in the pseudo-integrable billiards. 

If we consider an integrable system, A corresponding to an almost inte- 
grable system, P and let g be the genus of the surface R corresponding to P, 
denote by | A | and | P \ the respective areas of A and P, Gutkin [9^ proved 
that there exists a constant Ci such that 

F{l) = c,^^P + 0{l). (4.3) 

The constant Cie[l, M]. 
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4.2 Law of Proliferation of Periodic Orbits 



in Pseudo-integrable Billiards 

Our results in the previous chapter clearly showed that there are countable 
number of families of periodic orbits. By a family of periodic orbit, we mean 
an isolated trajectory closing after an odd number of reflections, or a band 
of trajectories closing after an even number of reflections. We have also seen 
that the periodic trajectories only occur in bands for the vr/S— rhombus bil- 
liard due to its equivalence with the barrier billiard. The number of families 
of periodic orbits of length less than x is finite for any x;we shall call this 
number by Counting Function, F{x) now on. 

In the subsections below, the law for an almost integrable and a pseudo- 
integrable billiard (TUl IH] is derived. 



4.2.1 The 7r/3-Rhombus Billiard 

We now present number-theoretic arguments to obtain F{x) for the vr/S- 
rhombus billiard. Subsequently, we shall discuss and compare our results 
with the above-mentioned results by Katok 3J and Gutkin The length 
of periodic orbits in the given family corresponding to the lattice point {q,p) 
where q,p are co-prime is given by 



l = C2L^{q^ + 3p^) (4.4) 

where C2 depends on family of periodic orbits as seen in the Table 3-1 of 
chapter (3). If / < x for a given family of points {q,p) should be counted 
in F{x). We can draw a circle of radius /i(= q^L) then, all points {q,p) 
having family of periodic orbit with length C2/1 < x (or li < xjc'i) should 
be considered for the calculation of F{x). Referring to Fig. 4.1, area of a 
quarter circle is (7rZ^/4) and the area of the square OABC is (/^). We shall 
denote the integer (fractional) part of a number by [...]({...}). The number 
of lattice points in OABC is 



= (go + l)Vv^-(go + l)^^ 



[ fa+D ] (4-5) 
On an average, [(go + l)/"\/3] is 1/2 (obviously); hence 
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iV^4 + ^^-^^ (4-6) 
Therefore, number of lattice points in a quarter circle is 

AT, = N{7^lllA)lll = ^ (4.7) 

Since the probability that two randomly chosen numbers are co-prime is 
(6/7r^) the number of co-prime lattice points is 



Nc= {-,){—) -^ + ^ (4.8) 



6 vrA^ V3ll Ay/3 
.^)(— ) ~ 27rL ^ 1^ 

where 

A = (4 - V3)/i/L + 2 - Vs. (4.9) 

For reasons discussed in the previous section, we are concerned in finding the 
number of points in a (vr/3) sector, the area of which is two-thirds that of a 
quarter circle. Hence, for (7r/3)sector, 

Taking only the dominant contribution (0(Z^)), with the help of Table 
(3.1) (cf. chapter 3), we can write for the number of periodic orbits whose 
length is < X as 

F{x') ~ 2(f oof.|,A^QQ|,(;.(2x) -|- Poocc-^oocci^-^^ ~l~ Poecc-^ oecciS^^ 

~l~-Poecc-^oecc(-^/2) ~l~ Peocc-^eocci^-^) ~l~ Peocc-^eocci^-^ / (^'H) 
~l~ -Pooce-^ooce (23^/3) -|- Pooce-^oocei^-^ / ~l~ Pooce-^ ooceiS^ / 

where [Poocc] is the probability that given co-prime lattice point is of 
odd-odd, centre-centre type and so on, NoocM) is total number of odd-odd, 
centre-centre type co-prime lattice points contained in the sector. Out of 
the four points, only one is odd-odd (or even-odd or odd-even), also two 
out of three points are of centre-edge type and one is of centre-centre type. 
Therefore, we have 



P = P = P =- 

^ oocc ^ oecc eocc „ 

y 
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2 

P = P = P = - 

^ ooce ^ oece ^ eoce g 

Thus we can write 



FW=2((4 + l + l + j + l + j)(5^) + (5 + - + 5)(5^)) 

or 

_ / 53 
~ V27y3vrL2 

In terms of |P| this can be written as 

F{x) = 0.049 733 (^^j x'^. (4.12) 

Of course, apart from the dominant term that is quadratic in x, there 
will be terms of 0{x) and 0(1). It is, however, important to note here that 
term of 0{x) is not related to the orbits periodic after an odd number of 
reflections ( hence, isolated). These terms only present more exact expression 
for F{x) arising from the above arguments. Their origin is in the points 
contributing to F{x) lying on the boundary of the sector. On the same lines 
as above, terms of 0{x) and 0(1) are found to be (26/8l7r)(4V3 — 3){x/L) 
and (12/277r)(2-\/3 — 3) respectively. Thus the counting function is of form 




4.2.2 The H-M Barrier Billiard: 

We now present our calculation for the Hannay-McCraw billiard Recall 
the discussion about this billiard in the previous chapter where we have 
argued that the periodic orbits in this system are the ones which hit any 
lattice point {q,p) in this array of barriers, the gradient of such trajectories 
will be given by \p/q\- Therefore we consider only the pairs(g,p) such that 
q and p are co-prime since they only give a primitive periodic orbits, and 
points {mq,mp), where m ^ Z, gives m repetitions of a primitive periodic 
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orbit corresponding to Each such {q,p) gives different number of bands 

or famihes of periodic orbits, depending on whether the pair is odd-odd (o, o) 
or even-odd (e, o) or odd-even (o, e). The length of periodic orbit in a given 
family corresponding to a lattice point (g, p) is given by 

/ = cL^{q'^+p^) 

where c depends on the number of families of periodic orbits. It can be seen 
from table (3.2) of chapter 3, that, for {q,p) as [12] 



4 one band, closing at (4g, 4p) for (o, o) 

2,2 two bands, each closing at {2q,2p) for (o, e) 

1,1,2 two bands closing at (g,p) and one at (2g, 2p) for (e, o) 

(4.14) 

If / < X for a given family of points (g,p),the contribution from this family 
of {q,p) should be counted in F{x). Drawing a quarter circle of radius h in 
the quadrant under consideration, all points (g, p) having family of periodic 
orbits with length cli < x (or h < x/c) must be considered for the calculation 
of F{x). The quarter circle is inscribed in a square OABC with side length li 
( Fig.4.2). thus the area of this square is If and the area of the quarter circle 
OAC is irll/A. We shall denote the integer (fractional) part of a number by 
[ ] ({})• The number of lattice points in the square OABC is given by 



N 




Taking fractional part of li/L, on an average as 1/2, we can write 

2 




Then the number of lattice points in a quarter circle is just 

(4.15, 



Since the probability that two randomly chosen numbers are co-prime is 
(6/7r^) the number of co-prime lattice points in quarter circle is 
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p 




Fig. 4.2: Lattice structure ofH-M billiard and single 
slit rectangal billiard is same though fundamental 
regions are not. End points of the barriers can be 
assigned integer lablelling. 
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{2^) ^1 + (ill) ^1 + If 



(4.16) 



The counting function can now be written explicitly as (using equa- 
tion (jHHl)) 



F{X) = PooNoo 



X 



+ 2PoeNo 



I) + 2PoeNoe {X) + PeoNeo (| 



(4.17) 



where again e.g., Poo is the probability that a given co-prime lattice point 
is of {odd, odd) type and Noo{x) is the total number of odd-odd co-prime 
lattice points contained in the quarter circle of radius x{= Nc). Trivially, 



Po. 



p = p 

-'- np. Pa 



1 
3 



thus 



F{x) = i [iV, (f ) + 3iV, (f ) + 2iV, 



V 327rL2 ) 



3_ 



(4.18) 



(4.19) 



This is the asymptotic law of proliferation of periodic orbits for system un- 
der consideration. How fast actual F[x) converges to equation ()4.19p depends 
on the rate of convergence of Poo, Poe, Peo and in accordance with equa- 
tions ()4.18j) and ()4.1fij) respectively. It can be easily seen that P converges 
rapidly to 1/3. In Table (4.1), we compare the actual number of co-prime 
pairs with the results obtained by equation ()4.16|) . It can be seen that (even) 
at X = 50, the % difference between the numerical and analytical results is 
only 2%. For similar reasons, we get an equally remarkable agreement in the 
case of the 7r/3-rhombus billiard as seen in Fig. (4.3). 



4.2.3 The Single Slit Rectangle Billiard 

This billiard has been already described in the last chapter. Again we con- 
sider only the pairs(g,p) such that q and p are co-prime since they only give 
a primitive periodic orbits, and points {mq,mp), where m & Z, gives m rep- 
etitions of a primitive periodic orbit corresponding to (g,p). Each such (g,p) 
gives different number of bands or families of periodic orbits, depending on 
whether the pair is odd-odd (0,0) or even-odd (e, o) or odd-even (o, e). The 
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Fig. 4.S:A log-log plot of the counting function, F(x) vs. the 

length I of the primitive periodic orbits. Dashed curve 
represents theoretical results eq. (4.13) and the solid curve 
represents actual F(x) obtained numerically. 
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length of periodic orbit in a given family corresponding to a lattice point 
{q,p) is given by 

where c depends on the number of families of periodic orbits, which is now 
changed due to a simple modification of H-M billiard to single slit rectangle, 
and is now given by (ref. Table 3.3 of chapter 3) 

{4 two bands, each closing at (4g, 4p) for (o, o) 
4 two bands, each closing at (4g, 4p) for (o, e) . (4.20) 
2,2 four bands each closing at (2g, 2p) for (e, o) 

Remaining analysis is same and will not be repeated here. Proliferation law of 
periodic orbits can be deduced similarly and will be of same form ax"^ + bx + c. 



The similarity of equations (j4.13p and (j4.19p indeed suggest an immedi- 
ate generalization. However, first of all, let us discuss the reason underlying 
the difference between the quadratic coefficient obtained by us and the one 
by Gutkin in the theorem paraphrased earlier in this section. The calcula- 
tions show that this difference is of an order or so. The reason is as follows. 
In considering the number of lattice points formed by stacking fundamental 
region of the corresponding integrable system, the condition of co-primality 
was not taken into account by Gutkin. As explained earlier an orbit labeled 
by a pair {mq,mp), where m & Z, gives m repetitions of a primitive peri- 
odic orbit corresponding to point {q,p) where q and p are co-prime. Hence, 
ignoring the co-primality condition leads to an over counting of the peri- 
odic orbits. Further, due to symmetry in the tessellated two dimensional 
plane, calculation need to be performed for (7r/4) sector in the single slit or 
H-M billiard and for (vr/G) sector in the 7r/3-rhombus billiard. In general, 
of course, for a domain with a discrete symmetry of order A^, only a (tt/N) 
sector needs consideration. Finally, one must note that a basic difference 
between the lattice generated by fundamental polygonal billiard and corre- 
sponding integrable system, lies in the incomplete tessellation of the plane 
by non-integrable billiards. For instance, the barriers are of zero width and 
finite length in the two examples considered in this section. It is this struc- 
ture that enables us to completely classify the orbits via integer labeling. 
The relative weight {Poo, Poe, Peo) in H-M billiard and {Poocc, Pooce) etc. in 
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the (7r/3)-rhombus billiard for different types of co-prime lattice points differ 
in different system and lead to a different quadratic coefficient. Hence, to 
give a general formula for the law of proliferation of periodic orbits exactly 
demands a complete enumeration and classification of periodic orbits. Al- 
though this important question cannot be answered today, a general recipe 
in the following is presented which comes very close to an exact formula for 
quadratic coefficient. 

4.3 Proliferation of Periodic Orbits Consid- 
ering Repetitions 

We now discuss the nature of this law if repetitions of the orbits are also 
counted, this finds application in the semi-classical theory of spectral statis- 
tics. If we follow on the similar lines as above, we obtain the asymptotic 
behaviour detailed below. 

If we are counting repetitions of the primitive orbits of length, I < x, 
then the primitive orbits of length between x and x/2 will not be repeated; 
the primitive orbits with length /; x/2 < I < x/3 will be repeated once; the 
primitive orbits with length x/3 < I < x/A will be repeated twice and so on. 
Thus taking into account these repetitions, one can write an expression for 
number of ' effective ' co-prime lattice points within a quarter circle of radius 
X, N{x) as. 



Nrix) = N,ix)-N, (I) +2 (iV, (I) - iV, {l))+-+n (iVc (^) - iV,(l)) 

(4.21) 

Here n is the largest integer less than x, we have neglected N{1){1 < 1) 
since there are no periodic orbits of length less than 1 in the system we have 
considered. Equation 1)4.211) can be rewritten as 



NJx) = NJx) + N, 



+ ... + N, 



nNJl] 



(4.22) 



Then, the modified counting function, Fr{x), for say H-M billiard becomes 



Frix) 



X 



Nri - ] +3Nr{ - ] + 2Nr(x 



X 



(4.23) 
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substituting equation (j4.22|) in equation (|4.23p . we get 

+ + ^ + (A) x-{n, + n, + n) (f )) 

(4.24) 

where n4 and n2 are the largest integers less than x/4 and x/2 respectively. 
A and B are given by 

Asymptotically {n — >■ cxd), ^(r;,) = 7r^/6 and = log(n) + 7 where (7) is 

the Euler-Mascheroni constant, equal to 0.5772157 The proliferation law 

with repetitions for other billiards considered above can similarly be deduced. 

The rate of convergence of actual Fr{x) to equation ()4.24|l depends upon 
the rate of convergence of actual Nr{x) to equation ()4.22p . In Table(4.1), we 
compare actual number of 'effective' co-prime pairs Nr{x) with that given 
by equation ()4.22p . Note that the % difference between actual 'effective' co- 
prime numbers and those obtained from equation ()4.22|l at a; = 50, is 5% ( it 
may be recalled that this figure is almost 2.5 times the one observed for Nc). 
It is for this reason that the convergence to the quadratic law (equation ()4.24|) ) 
is much slower if one considers repetitions. 

To conclude, we have shown analytically that the rate of proliferation of 
the periodic orbits is exactly quadratic (in length) asymptotically. We have 
seen that this is in complete agreement with extensive numerical calculations 
on two model pseudo-integrable systems. The reason underlying the asymp- 
tote, ax'^ + bx + c, to the counting function is clearly related to the tessellation 
of the two dimensional plane by the fundamental region of the billiards. 

4.4 Generalization 

It is well known that a rational polygon can periodically tile a surface ev- 
erywhere flat, in the sense of null Gaussian curvature, except at isolated 
vertex points of singular negative curvature. A periodic structure that tiles 
the almost everywhere fiat surface may consist of several polygons and hence 
space can be assigned distinct labels (albeit complicated) taking account of 
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different periodicities, in a spirit similar to the one presented above. To enu- 
merate distinct primitive periodic orbits, one needs a condition analogous 
to the co-primality condition required by the two systems discussed above, 
since out of all lattice points lying on the same line of a given slope only one 
will give a primitive periodic orbit. Let us denote the probability of the "co- 
primality condition" to be satisfied by distinct labels Pc- Furthermore, the 
classification entailing each distinct label will give rise to relative weights in 
which the orbits will be distributed, let us denote it by Pj{j denote classes). 
For a polygon with symmetry group of order iV, the points to be considered 
will be restricted to an/N sector. This number can be written as 

i 

where the summation is over all the periodicities. And number of "co-prime" 
points are Nlc{x) = PcNl{x). For each class of periodic orbit for which the 
weight is Pj, there may be k types of periodic orbits closing at length ^kjX. 
We assume that k remains constant for given class. With this the counting 
function becomes 

j k V'S'JJ/ 

and hence coefficient of F{x) — ax^ + bx + c are 

a = ^E^iEPjE^ik 
b = ^i:P^EP,E^ik 

I j k 

c = ^EiiEPj 

Hence, counting function will quadratic and convergence will depend on 
both Pc and Pl- However, If k is not constant but depends on x then above 
analysis will have to be modified and F{x) may not be quadratic. 
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X 


Nc 


Nc 


Idiff 


Nr 


Nr 


/dtff 






e.q.\A.\{\\ 


Col.2,3 






Col.5,Q 


3.0 


2 


5 


150.0 


2 


4 


100.0 


4.0 


6 


9 


50.0 


6 


9 


50.0 


5.0 


10 


14 


40.0 


12 


17 


41.6 


6.0 


16 


19 


18.7 


18 


23 


27.7 


7.0 


20 


26 


30.0 


26 


34 


30.7 


8.0 


28 


34 


21.4 


36 


46 


27.7 


9.0 


36 


42 


16.6 


48 


59 


22.9 


10.0 


46 


52 


13.0 


62 


75 


20.9 


20.0 


190 


200 


5.2 


280 


313 


11.8 


30.0 


424 


443 


4.4 


654 


707 


8.1 


40.0 


764 


782 


2.3 


1188 


1264 


6.3 


50.0 


1192 


1217 


2.0 


1876 


1977 


5.3 


60.0 


1718 


1747 


1.6 


2720 


2848 


4.7 


70.0 


2330 


2372 


1.8 


3724 


3876 


4.1 


80.0 


3048 


3093 


1.4 


4884 


5063 


3.7 


90.0 


3876 


3910 


.8 


6208 


6406 


3.2 


100.0 


4770 


4822 


1.1 


7684 


7912 


2.9 


200.0 


19088 


19193 


.55 


31066 


31593 


1.7 


300.0 


42972 


43114 


.33 


70162 


71009 


1.2 


400.0 


76402 


76585 


.24 


124974 


126145 


.94 


500.0 


119372 


119604 


.19 


195484 


197005 


.78 



Table 4.1: Comparision between actual and theoretical values for and Nr 
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Chapter 5 



Two-Point Cluster Function 
and Form Factor: A Diagonal 
Approximation 

5.1 Introduction 

We Consider the discrete spectrum Ei, E2, E^, ■ ■ • of a /-dimensional bound 
quantum system. As mentioned in the chapter (1) one can define a "global" 
function N{E) which gives number of levels less than or equal to E on the 
set of this discrete levels. The function N{E) is obviously a staircase like 
function, since it jumps by unity at the occurrence of an eigenvalue Ei and 
remains constant otherwise. Of course if spectrum consist of degeneracies 
then jump in N{E) at the particular degeneracy will not be unity but equal 
to number of degenerate levels at particular point on energy axis. Typical 
functional nature of N[E) is therefore non-analytical. When we approximate 
N[E) by a semi-classical expression we use approximate analytical form for 
the exact staircase. We use same symbol N{E) for both actual staircase and 
its semiclassical approximation, mostly we will be talking about semiclassical 
approximation below. 

As mentioned earlier it is possible to separate semi-classical approxima- 
tion to N{E) in a smooth part Nav{E) and the fluctuating part Nosc{E) (i.e. 
N{E) = Nav{E) + Nosc{E)). In the particular case of a particle (of mass 
m — 1/2) in polygonal box (or polygonal billiard) of area A^, Nav{E) is 
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given by 

Na.iE) = ^-^ + C (5.1) 

where is length of the parameter of the boundary, and C is a constant 
containing complex information on the geometrical and topological properties 
of the billiard. 

This average behaviour may be eliminated in order to characterize and 
compare the fluctuation patterns of different systems whose corresponding 
Nav{E) is not same. For this purpose it is convenient to ''unfold'' the 
original spectrum {Ei} through the mapping E i — > r 

n = Nav{Ei). (5.2) 

The effect of fl5.2|) is that the sequence {rj} has on the average a constant 
mean spacing equal to unity, irrespective of the particular form of the function 
Nav{E). By construction rj^z — l/2(z = l,2,...) and the departures 

5i = r,-{i- 1/2) (5.3) 

of Tj from its average value i — 1/2 are the level fluctuations. As a result 
of unfolding the spectrum in a way mentioned above it is obvious that the 
average part of level density of unfolded spectrum dav{y) is unity. 



5.2 Spectral Fluctuation Measures 

The level fluctuations can be characterized in a systematic way using the 
A;-point correlation functions and measures derived from them. The fc-point 
correlation function is defined as |2] 

I oo oo 

Rk{Xl, ...,Xk) = , J ... J PniXi, ...Xn)dXk+l...dXn. (5.4) 

— oo — oo 

where P„ is a joint probability density function of the levels {xj}j=i...„ and 
1 < k < n. The Rk is thus a probability density of finding a level (regardless 
of labeling) around each of the points Xi...Xk, the positions of the remaining 
levels being unobserved. Each function Rj.; k > 1 contains terms of various 
kinds describing the grouping of k levels into various subgroups or clusters. 



94 



In practice it is convenient to introduce fc-level cluster functions obtained 
from by subtracting out the lower-order correlation terms, as 



T,(xi...Xfc) = ^(-l)'=-"^(m-l)!n^G,(xt,with k in G,). (5.5) 

G 3=1 

Here G stands for any division of the indices [1, fc] into subgroups [Gi, Gm]- 
For example, 

Ti(x) = Ri{x) (5.6) 
and 

T2{Xi,X2) = -R2{Xi,X2) + Rl{Xi)Ri{x2) 

Measuring the energies in the units of the mean level spacing A = 1 and intro- 
ducing the variables yi = Xi/X, (this is nothing but the unfolding mentioned 
earlier) the cluster functions can be written as 

Yk{yi,---yk) = lim A"T„(a;i, ...Xfc). (5.7) 

The Yfc's are then well defined and finite everywhere. The cluster functions 
being isolated from the effects of the lower correlations, vanishes as the sepa- 
ration 

r =1 ri — r2 | becomes larger and larger. Among these two- level cluster 
function Y2 is of prime importance. When the collection of levels are treated 
as a classical Coulomb gas, Y2 defines the shape of the neutralizing charge 
cloud induced by each particle around itself j2]. As mentioned earlier many 
fluctuation measures can now be expressed in terms of the Yk 's. We shall 
discuss this at appropriate place. 

Consider now one point measure, the energy level density function 
d{y) = J2^{y ~ y-i)i where 5 is usual Kronecker's delta function. The aver- 

i 

age with respect to the probability distribution function P(?/i. ..?/„) (so called 
ensemble average) of the level density function is given by 

00 00 

{{d{y))) = J dyi-- - J dynd{y)P{yi...yn) 

—00 —00 

00 00 
= n J dyi - ■ ■ J dyj{y - yi)P{yi...yn) ' 

—00 —00 

= R,{y) = T,{y) 



(5.8) 
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Also the density-density correlation function (for energy) is defined as 

n n oo oo 

{{d{y)d{y'))) = E E I dyi ■ ■ ■ J dyj{y - yi)5{y' - yj)P{yi...yn) 

i=lj=l— oo — oo 

n oo oo 

= nj^ J dyi--- J dyj{y - yi)6{y' - yi)P{yi...yn) 

i=l— oo — oo 

n oo oo 

+n(n-l) E I dyi-- - J dyn5{y - yi)5{y' - y2)P{yi.. 

j^j=l— oo — oo 

= 5iy-y')R,iy) + R2{y,y') 
Hence, from equation (j5.6|) . we get 

my)d{y'))) = 6{y - y')R,{y) + Ri{y)Ri{y') - Uy, y'). (5.9) 

Since the ensemble average is same as the spectral average(see e.g.|3jj), which 
is defined as (for level density) 



Idiy)) = - dx d(x), (5.10) 
n J 

we can write in an asymptotic limit 

(%)) = Ri{y) = Y,{y) = daM = 1 (5.11) 

and 

{d{y)diy')) = 5{y -y') + l- Y^iy, y'). (5.12) 

Thus two-point correlation function R2{yi,y2) or cluster function 1^2(1/1,7/2) 
are related to the spectral average of density-density correlation function as 

R2{y,y') = {d{y)d{y'))-6{y-y') (5.13) 

and 

Y^iy, y') = l + 5{y - y') - {d{y)d{y')). (5.14) 

Equations ()5.13|1 and ()5.14|1 are the important relation to develop semi- 
classical understanding of spectral fiuctuation properties in which we are 
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interested. Since we can write density of states as sum of average and fluctu- 
ating part {d{y) = dav{y) + dosdy)) the density-density correlation function 
becomes 

{d{y)d{y')) = {{daM + dosc{y)){da,{y') + dosc{y'))) 

= dav{y)dav{y') + dav{y){dosc{y')) 
+dosc{y){dav{y')) + {dosc{y)dosciy')), 

since (dosdy)) = {dosdy')) - 0> and also dav{y) = 1, we finally get 

{d{y)d{y')) = 1 + {dosdy)dosdy'))- (5.15) 
Using this, equations ()5.13|) and ()5.14|) can be rewritten as 

R2{y, y') = {dosdy)dosdy')) + ^ - ^{y - y') (5-i6) 

and 

Y2iy,y') = 5{y - y') - {d^y) dosdy'))- (5-17) 
These are the equations that we will use below to study spectral fluctuation 
properties as we can substitute semi-classical expression obtained in chapter 
(1) for dosdy)- Note, existence of S{y — y') in eqs. ()5.16p and ()5.17p . Since 
i?2 or Y2 must be free of ^-function, density-density correlation function also 
must contain similar ^-function. Indeed this is due to the self-correlation of 
a level. 

All the relevant information to evaluate semi-classical periodic orbit sum 
is obtained in the previous two chapters. For an easy reference we will repro- 
duce briefly those results in the following section. Further the example, we 
consider here (vr/S-rhombus billiard) also possesses a symmetry with respect 
to both the diagonals. This results in modification of the semi-classical ex- 
pression for the density of states, since one can now work with decomposed 
state space due to symmetry. These issues are also discussed in the following 
section before we take up study of statistical measures. 

5.3 Semi-classical Density of States 

As stated in chapter (1) the oscillatory part of density of states in case of 
pseudo-integrable is given by|4| 

d{E) = -^P^ReY, J J dciH^Q^ {^^2mE^ exp(mi7r). (5.18) 
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where all symbols are as explained in chapter (1). Here we take m = 1/2, 
then by taking real part and in the asymptotic limit {h 0) we get (note 
that a's are integers) 

dosciE) = (^) E~y^ E ^ cos (|v^ + nia, - 1)) . (5.19) 

We now ' 'unfold' ' the energy using only first term of A^^^ for the rea- 
sons of simplicity (as a trade off between simplicity and accuracy), through 
mapping E ^ tq : 

- - si- 

This means unfolded levels will not have mean spacing equal to unity ex- 
actly but only approximately. In fact from equation (j5.1|) it is easy to 
see that mean level spacing after unfolding is less then unity by a factor 
(1 — Ld/ (^/Ai{A^r\i) + C/tq). Since for semi-classical analysis is valid for 
large value of tq, it is obvious that unfolding levels using just a leading term 
of Nav{E) will not introduce serious error in the analysis for large value of 

From equation (j5.2(J|) we can see that 

— = — ^ ~ dav{E), 



dE Anh' 



and hence, 



d{ro) = d{E) (^^^ c^l + dUE)[j§^ =l + dUro). (5.21) 

The oscillatory part of the density of states feq. ()5.19|) ) can then be written 
in terms of rescaled energies as 





dosciro) = -rh-ui X ^ cos + («. - 7)^ • (5-22) 



As mentioned above when a quantum dynamical system possesses a sym- 
metry, its state space may be decomposed into subspaces of definite symmetry 
type. The Schrodinger equation can then be restricted to these subspaces 
resulting in a symmetry- projected spectrum. It is interesting to know how 
these symmetries influence the fluctuation properties of the whole spectrum. 
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5.3.1 Symmetry-projected Green's function 

Let H he a. quantum Hamiltonian for a system of / degrees of freedom, 
invariant under a discrete group G. We will assume that G acts classically as 
a group of Euclidean point transformations ( i.e. combination of translations 
rotations, and reflections). The quantum mechanical action of G in the x 
representation is given by 

Uig)\^) = \g-^) (5.23) 

where g is an element of G and g ■ x denotes its action on x. By virtue 
of this symmetry, H may be restricted to a subspace which is invariant un- 
der the symmetry and whose states transform according to an irreducible 
representation of G. The projection onto an invariant subspace is given by 



"^""T.xMu^g), (5.24) 



\G 



where, Xm{g) is the character of the m*^ irreducible representation of dimen- 
sion dm, and \G\ is the order of group GjS]- 

A semi-classical approximation for Gm{E), the trace of the symmetry- 
projected Green's function is defined as 

GUE) = TT[Pm{E-H)-']. (5.25) 

The poles of gm{E) are the energy eigenvalues of symmetry m and the pole 
strength give the eigenvalue degeneracies. One can now follow steps given 
in chapter (1) to obtain the periodic orbit sum formula. The symmetry- 
projected Green's function, Gmi^,^'] E) = {x\Pm{E — H)~^\x.') can be ob- 
tained by restricting to an invariant subspace only which is given by^ 

G„(x, x'; E) = ^J2 xMGig ■ X, x'; E). (5.26) 

l^^l g&G 

It may be noted that the full semi-classical Green's function involves sum over 
orbits between x' and x, its symmetry projection involves a larger family of 
orbits. The sum over classical orbits is taken for orbits that begin at x' and 
end at (7 ■ X, a point related to x by symmetry. These orbits can be associated 
in a one-to-one fashion with orbits on a symmetry reduced phase space. 

A general form of symmetry (discrete) projected oscillatory part of density 
of states can be obtained from the corresponding trace formulajU]. In the 
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case of pseudo-integrable polygonal billiards it reduces to (in an unfolded 
spectrum) 

JJ^. j E ^ COS — - - + /.,vr j xmign (5-27) 

Ft J j 

Here, the terms Ip, Ap are length and band area of the periodic orbits in the 
primitive polygon that tessellate full phase space. The sum over all m, of 
equation ()5.27|) reproduces equation ()5.22|) . 

The spectrum of the vr/S-rhombus billiard is composed of superposition of 
four different modes due to system symmetries j7j. Semi-classically, it is not 
possible to distinguish all four modes Q separately. However, it is possible to 
separate equilateral triangle modes(odd parity classes) where wave function 
is zero on the shorter diagonal and pure rhombus modes(even parity classes) 
where gradient of wave function is zero on the shorter diagonal. This can 
be done by incorporating and exploiting the symmetry of the system in a 
semi-classical treatment, to obtain symmetry projected spectrum. 

We will follow treatment of [6^ to study effect of symmetries on the fluctu- 
ation properties of the spectrum. We take longer diagonal along y-axis and 
shorter diagonal along x-axis. 7r/3-rhombus billiard has symmetry group 
C2<S>C2, consisting of reflections about x and y axes. The elements of the 
group can be denoted by x°-y^, where x and y represent commuting reflec- 
tions about a respective axis. The primitive cell ^ jH| is then half equilateral 
triangle having sides half of longer (Cy,) and shorter (C^) diagonals and one 
side of the rhombus C. The group has four one dimensional {dm in ()5.27|) 
is therefore 1) irreducible representations, labeled by numbers {p,q). These 
numbers describe parities under reflections about x and y axes respectively. 
For states, even under x reflection, p = 0, and odd under x reflection p = 1. 
Similarly, q takes values or 1 depending on parity under y reflection. The 
characters of the group can be easily determined in this case and are given 
as 

XvM) = exp{-7™(pai + qhi)} (5.28) 

where, a/, hi are number of times Z**^ primitive periodic orbit bounces from £3, 
and Cy respectively, n counts repetitions of periodic orbits and gf in group 
element x"'?/^'. Pure rhombus modes are therefore represented by p = 0; g = 
or 1 and the equilateral triangle modes byp=l; g = Oorl. 

Recalling the comments about the Maslov indices in the chapter (1), one 
can see that the character of the group actually play a role of the Maslov 
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indices n in this case. Thus, for the pure rhombus modes Neumann boundary 
condition is apphed to the shorter diagonal and then bounces on the same 
are not taken into the consideration. For equilateral triangle modes however, 
since one have Dirichlet boundary condition on Cx, bounces on it should be 
counted. 

From the polar constructions we have used in chapter (3) to determine 
and classify the periodic orbits a/ and bi can be determined easily. It can be 
seen that bi is always even, hence qbi does not play any role here. This is the 
basic reason that one can not separate these two parity classes (belonging to 
integrable and pseudo-integrable parts) further into the four different modes. 
Also it turns out that for all the centre-edge orbits, ai is even. However, 
for all the centre-centre orbits there are three bands of periodic orbits each 
occupying area of AAr {Ar is area of V). Two of these bands bounces odd 
number of times on Cx and one band bounces even number of times (i.e. ai 
is odd for two bands and even for one band). The same is true for side £, 
since total number of bounces on three sides of V are always even. Using all 
these facts discussed above, equation (j5.27j) becomes 

'^-'^^'■^ = ( j/4^3/4 ) E ^ COS - ^ + f^j^ - (5.29) 

Note, repetition index n is absorbed in aj here. Aj,lj and Ar represent 
respectively the band area of periodic orbits of primitive cell V, length of 
the periodic orbits and area of half equilateral triangle. The Maslov indices, 
IJ,j here counts bounces only on the side C of the primitive cell V. Thus for 
the pure rhombus modes since p = 0, fij will play important role, while for 
equilateral triangle modes (p = 1), (fij — aj) is always even for all the periodic 
orbits. If therefore, /I's are neglected one will end up in erroneous results. 

Using above discussion and following the steps of IE] we write density of 
states for the pure rhombus modes as 

for centre-edge orbits. For all centre-edge points, Aj will be 12Ar i.e. same 



101 



as phase space area. And for centre-centre orbits we write 




(5.31) 

+ ^'j jh cos (/, - f ) - E ■ jh cos (/,V5 - f ) 

3 j 

Here, superscripts a, e, o over summation sign represents all repetitions, 
even repetitions and odd repetitions of the primitive periodic orbits. This is 
simple to understand, since centre-edge points do not result in periodic orbits 
of different lengths in case of the rhombus there contribution will be same 
as the rhombus or in the half equilateral triangle. However, note that each 
centre-centre point results in two kinds of bands of the periodic orbits one 
having the same length as that of the periodic orbits in the half equilateral 
triangle and other having length twice of the same. These orbits do not 
belong to the half equilateral triangle. The density of states for equilateral 
triangle mode has same form as ()5.30p . It is simple to check that adding 
these contributions we do get density of states for the complete rhombus. 
We shall now turn our attention to study of two-level cluster function. 

5.4 Two-Level Cluster Function 

In this section we shall obtain a closed analytical expression for two level 
cluster function in pseudo-integrable billiards using semi-classical framework. 
To do this it is clear from equation (|5.17|) that we need to evaluate semi- 
classical product {dosc{y)dosci3j')) , (i-e. spectral average of density-density 
correlation function. 

5.4.1 Density-density correlation function 

As defined above density-density correlation function (oscillatory part) is 
given by dosc{y)dosc{y')- We consider correlation of levels around tq i.e. we 
will put y = To + ri and y' = + r2 : ri, r2 « then using equation (j5.22|) 
correlation function becomes 

— 1/2 

dosc{ro + ri)dosc{ro + r2) = ^"^^3/2 J2j,k jrr^i {cos (^jV^o + ri + tt (a^ - i))) 

X cos (y^ (/fc V^o + r2) + tt (a^ - i))} 
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or by using simple trigonometrical identities 



dosc{ro + ri)d 





{sin {IjV^o + ri + hV^V+r^) + '^^ 

+ cos {IjV^o + ri - hy/ro + rg) + } 

(5.32) 

where k± = n^aj ± a^), terms due to Maslov indices(i.e. number of bounces 
on the billiard wall). In the case of polygonal billiards it can be seen that 
the number of bounces of a particle on wall is always even if boundary con- 
ditions on the all segments of the boundary are same either Dirichlet or 
Neumannjn]. Hence, we will drop them from our analysis until we consider 
an example where they becomes important (these are the cases where bound- 
ary conditions are mixed i.e. on some segment it is Dirichlet, on the others 
it is Neumann). 

It is easy to see that due to local averaging the first term in the above 
equation will be negligibly small even when statistically significant number 
of levels are considered. The second term also vanishes, unless 



< 1 



The reason for this is that in the semi-classical limit this term will oscillates 
rapidly fOlE] (it is well known that if argument of cos or sin function is large, 
any small change in it will result in rapid oscillations). Due to these rapid 
oscillations, spectral averaging will result in negligibly small contributions 
unless above condition is satisfied. 

Since r << ro we can approximate i/ro ± r ~ v^(l ±r/2ro) to write 
above condition as 

' ' r 



R 



5 + — / 
2ro 



< 1 



(5.33) 



where for a simplicity we assume with out loss of any generality Ij = I + 6/2 
Ik = 1 — 6/2, and ri(2) = +(— )r/2. Note the presence of correlation range 
'r = |ri — ' in the denominator of the second term, which appears because 
we are dealing with density-density correlation function directly and not with 
it's Fourier transform as generally done. The equation ()5.33|) indicates that 
one can separate the correlation function ()5.32j) in two parts as 



dosciro + ri)dosc{ro + rs) = J2 



U—lh 



E 
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The first part is known as diagonal part {S = 0) and other part is called as 
off-diagonal part (6^0). 

Now let us take close look at condition ()5.33|) which can be written as 



< ^ 

This is equivalent to following conditions 



r 

5 + — / 
2ro 



A 



R 







for0<5: 0<5<y^-4; (5.34) 



for5<0and |5| < ^ : ^1-\^<\5\<^1 



for 5<0 and |5| > ^ : ^/ < 1^1 < + 

One can see that in lim tq —>■ oo, (or equivalently h —>■ 0),: 6 —>■ 0, hence, 
off-diagonal contributions may not be important in the asymptotic limit and 
one can safely rely on the diagonal approximation only. Generally, one is 
interested in the asymptotic properties of spectral fluctuations only. 

In practice, however, one always works with finite ro, in fact most of the 
numerical studies deals with at most few thousands of levels. It should be em- 
phasized that in such cases off-diagonal contributions may not be negligible 
since finite number of levels are considered. 

When 6 = (i.e. Ij = 1^) this condition becomes 

' (5.35) 




2'Kr 



where Ih = a/^ttArTo is ' 'Heisenberg length' ' (which corresponds to well 
known Heisenberg time tn = 27ihdav{E) ,via relation 1^ = 2tj{VE)- 

The equation (j5.35p stipulates a condition for the validity of the diagonal 
approximation. For any given Ih, it is sufficient to work with diagonal terms 
for the length of periodic orbits satisfying above condition, since then one 
can see that in general 5 ~ for sufficiently large number of levels. This 
condition is similar to the one obtained for diagonal approximation for Fourier 
transform of Y2 (known as the form factor), except the presence of correlation 
range r in the denominator. It is obvious from this condition that for any 
fixed Ih diagonal approximation is strictly valid only at r ~ 0. For << r 
and In/'^f' < I the contribution of off-diagonal terms may not be insignificant. 

Here, in this and next chapter we consider diagonal terms only. Our 
results will therefore valid for r << 1. 
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5.4.2 Diagonal approximation 

We can neglect first term of eq. ()5.32j) . Since we are interested in correlations 
over the number of few mean level spacings, much smaller than total number 
of levels considered (r << ro) we can approximate y/r^ ± r by yro(H-r/2ro) , 
then density-density correlation function becomes 



where r = |ri — and Ih = \/47rA/jro, the Heisenberg length. 

We have seen in chapter (3) that the periodic orbits occurring in the 
pseudo-integrable polygonal bilhards can be classified into different families 
according to the nature of lattice points corresponding to the periodic orbit as 
well as to the projected phase space area occupied by the bands of periodic 
orbits. We have seen this via few examples of such billiards. For these 
billiards density-density correlation function can be written as 



Here, subscripts denotes classes of bands of the periodic orbits, 
are projective phase space area occupied by the periodic orbits in the given 
class, Qa 9(3 denotes degeneracy in the lengths of the periodic orbits. Other 
symbols have their usual meaning. The second term on RHS takes care of 
systematic length degeneracy occurring in different classes. However, we do 
not take into account accidental degeneracy here. The summation over j in 
the above equation can be replaced by integration in the continuum limit, 
choosing a proper measure. Heuristic arguments to obtain such measure are 
given below 

5.4.3 Summing over all the periodic orbits 

Consider a length interval (/, / + dl) such that dl « I but sufficiently large 
enough to contain many events (by events here, we mean occurrence of a 




(5.36) 





(5.37) 
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periodic orbit of a length within the above interval). The proliferation law 
of the periodic orbits, which gives the average number of periodic orbits of 
length < / for a class a, has a form 

Fail) = aj^ + Oil) 

where aa is system dependent constant. In our analysis we will neglect 0(Z) 
term. Number of periodic orbits having length between (dropping subscript 
a temporarily), I and / + dl, is given as n=2aldl. Hence the mean spacing 
between periodic orbits is l/2al. One can divide interval dl in n cells each of 
the length equal to the mean spacing. The length of the i*^^ periodic orbit, 
he (/,/ + dl) can be written as li=l + i/2al + r]i. Here, r^j's represent local 
fluctuations in the actual periodic orbit lengths from the average lengths 
/ + i/2al. It is therefore, reasonable to assume that ?7i's are distributed 
symmetrically around zero. Most of the rji' will be closely distributed around 
zero with small fractions extending up to few mean level spacings. The 
summation of form J2 cos(^Z)// over j for Ije (Z, I + dl) can be written as 

g cos(^/,) _ g COS [Ail + ^ + r^j)) / J 

j=o h j=o ^ ^ ^'^^^ ^ 



In our case •^^^Ja^''^ ~ '^'^'^I^h << 1, hence Arjj ~ 0. One can also notice 
that 2nd and 3rd terms in the RHS bracket of above equation are << 1. 
Using these facts one can write leading order terms of the above summation 
as 



!^i cos [A [l + ia]) cosiAl) / j \ sin(^/) . / j 

> ^ ^ — ~ ; > COS 7 ; > Sm 

f^o h I J^o ^2a/y / fr'o \2al 

Now summations on RHS can be easily performed by substituting, n 
2a I dl and using trigonometric relations (see e.g. [12] page 30) 

V^n-l ^°^iM+^)) ^ cosjAl) ^^^.( A\Adl I sin(^/) Adl 
^j=0 — I COLI^^^^J 2^1 2 

^ cos(^;) fry Y^oo / -I \k 2^''B2k A I sin(AO Adl 

— J [zaiaL) 2^fc=ol '-J (2ky. Aai ~^ i 2 



^(2a/d/) \l-Sr^+Oif) 



sinjAl) Adl 
I 2 
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Here, B2kS are Bernoulli's numbers. Neglecting terms of order 0{l ^) and 
higher as their contribution to summation is negligibly small, we get 



n-l 



COS 



^ -(■^('+^)) ^ £^2a/d/ (5,38) 
Same measure can be deduced from following simple relation 



However it should be noted that 2nd term will have non-trivial contribution 
unless ^ << 1 as indicated above. This is the reason for which one can not 
apply similar treatment to periodic orbit sum in density of states(e.g. eq. 
()5.22|) ) or for off-diagonal terms as well. 



Thus, returning to discussion on correlation function, using above results 
the leading term in density- density correlation ()5 . 3 7|) becomes 



dosciro + ri)dosciro + rs) 



1/2 



2 .3/2 {J2a Oa^a^al, 



a/3 



(5.39) 



where 



la — I dl COS — 

' ^min,a \ ^ U 



(5.40) 



I^y can be similarly defined and a^i, is appropriately chosen coefficient. In 
the subsequent analysis this second term will not be shown explicitly, though 
it will be included in the calculation of a specific example. The choice of a 
measure enables us to write two-point density correlation as 



dosc{ro + ri)dosc{ro + ^^2) = XI 



sm 



2-1^1 



max, a 



Ih 



sm — 



vrr 



vrr 



where 



R 



(5.41) 



(5.42) 
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and Ih has usual meaning as above. In limit /max ~^ oo first term in the 
bracket represents contribution to the self-correlation term (i.e., 5{ri — r2)) 
of the levels. In writing the two-level cluster function Y2 we remove this 
self-correlation singularity as in eq. ()5.19p . It may be noted that eq. ()5.4ip 
indicates that F2(ri,r2) is just a function of r = |ri — r2|. Hence, we write 
the cluster function as 

Y2{r) = Y: ^^a^ ^ (5.43) 

Here, (...) denotes spectral averaging as stated earlier. 

5.4.4 Taking spectral average 

The spectral average refers to all levels in an energy range that is classically 
small, i.e. small in comparison with tq, but large in comparison with the 
mean level spacing, which is unity in an unfolded spectrum or to be more 
precise large in comparison with the energy range defined by outer energy 
scale defined by Berry [lOj. The outer energy scale is given by h/T^i^, where 
Tmin is the period of the shortest classical closed orbit. It may be noted that 
each classical orbit with period T causes deviation(±) of a level from the 
mean position by order h/T. Hence, outer energy scale corresponds to the 
largest deviation of a level from the mean. In terms of number of mean level 
spacing it is given by Tout = ^dav/T^m ~ ////Svr/niin- In our case this range 
(say, a) then roughly becomes y/r^ << a << tq, since Ih ^ O (v^) • 
We now consider spectral average of a function / {At / ^/roj 

— CTjA 

Since a' << tq we can approximate this as 




Expanding / in the Taylor series as 
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If Ar << y/r^ we can approximate spectral average as 



Which can be easily evaluated and we get 

It is thus obvious that second and higher order terms can be neglected (since, 
a « tq) if /' and higher derivatives are also small. We will encounter here 
trigonometric functions (e.g. sin, cos and sin integrals) only. 



Thus after spectral averaging the two-level cluster function becomes 



sm ' • 



Y2{r) ^ ^ /""^ ^ (5.44) 

Trr 



As per convention negative sign of Y2{r) represents positive correlation 
and vice versa. For the finite spectrum occurring between rmax and rmin 
two- level cluster function is given by[Tni E] 

>2(r) = (5(r) - — / dro {dosc{rQ + n)dosc{ro + r2)) (5.45) 

''"min 

where Ar = r^ax — ''^min- Note, that we have used here a simple box function 
as a window function defined in J3] which takes care of the finite number of 
levels and should not be misunderstood as double averaging. Using eq. ()5.44|) . 
Y2{r) becomes 

Y2{r)^Y.^^A^ fdro sin f ^^^"l (5.46) 
a vrrAr 7 \ Ih J 

' mm 

Integration in ()5.46|) is thus easy to carry out and Y2{r) can then be 
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written in simple analytical form as 



a 



- r. 





+ cos ( 



cos 





r 




(5.47) 



where, recall that is Alglaa/irAR, Z„,max(min) is ^^^Ar rmax(min) and si(x) 
is sine integral — dt sm{t)/t. 

Apart from dependence on the energy window (rmin, Tmax) selected, one 
can see that Y2{r) depends on two parameters (1) /mm,a, shortest periodic 
orbits in each class a (we will call it length parameters) and {2)K.a, which is 
a measure of phase space area occupied by bands of periodic orbits in class 
a. Without referring to any system, it is interesting to study how these pa- 
rameters affect the behaviour of Y2{r). Fig. 5.1 shows effect of variation in 
length parameters where we have considered three classes of periodic orbits 
with same value of parameter k and energy window. The length parameters 
are, for curve a) 10,20,30, b) 10,20,40, c) 10,20,50, d) 10,20,70. The k's are 
same for all classes and curves: .3333 .3333 .3333. Number of levels consid- 
ered are 100-1000. Curve (g) is for G.O.E. As one increases highest length 
parameter one gets more and more level repulsion. This indicates that the 
presence of relatively longer periodic orbits will result in level repulsion in 
the quantum spectrum. In Fig. 5.2 we show effect of variation of parameters 
K for same values of length parameters. The k parameters are, for curve a) 
.6,.2,.2, b).2,.6,.2, c).2,.2,.6, d).l,.l,.8. The "lengths" are same for aU classes 
and curves: 10,20,30. Number of levels considered are 100-1000. One can see 
here that as weightage k increases for the longer length parameters, Y2{r) 
deviates from Poissonian (i.e. Y2{r) =0) showing more and more negative 
correlations or level repulsion. As far as the effect of energy window is con- 
cerned we can see from above equation that as r^ax — ^ oo and r^am <^ '"max 
so that Ar ~ Tmax, Y2{r) ~ maXa(47r/inin,a/ V^max), hence approaches a Pois- 
sonian as Tmax ^ oo in the diagonal approximation. Rate of approach 
is to Poissonian is basically governed by l/^/r^^. It may be noted that 
maxQ,(Zinin,a) IS also important parameter since it will decide when one would 
expect Poisson like results. For example in the case of 7r/3-rhombus billiard 
(pure rhombus modes ) meiXa{lrD.m,a) — 2\/93C, {C being side length of the 
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Figure 5.1: Two level cluster function: effect of variation in length parameters 
(see text) 
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Figure 5.2: Two level cluster function: effect of variation of parameters k 
(see text) 
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billiard), hence to see real asymptotic behaviour one should consider num- 
ber of levels must be at least a order of magnitude larger than 1600 levels. 
Unfortunately there exist few numerical results in this range. Most of the 
numerical experiments can therefore be termed as sub-asymptotic. In such 
cases as mentioned above off-diagonal contributions may be significant. Any 
comparison of our results with numerical experiments should be made by 
keeping this discussion in mind. 

We now consider two-point cluster function for vr/S-rhombus billiard. We 
reproduce relevant information in table 5.1 and 5.2. In Fig. 5.3 we show 
Y2 for complete 7r/3-rhombus bilhard, including both even and odd states 
that corresponds to two pure rhombus modes and two equilateral triangle 
modes. In Fig 5.4 we consider two pure rhombus modes combinedly (i.e. even 
states). Effect of number of levels considered is shown via curves a)represents 
energy window (47,370) b)(420,740) c)(47,743) and d)(1000, 2000). It is 
evident from these figures that even states are more ' 'away' ' from the 
Poissonian than the complete rhombus case. This is obvious due to fact 
that complete rhombus billiard spectrum is superposition of spectrums of 
four modes described above. One can also see that as we increase number 
of levels Y2 approach Poissonian. Since spectral measures we consider here 
directly depend on I2 we expect similar trend in case of these measures. 



5.5 The Form Factor 

In this section we will consider one of the interesting quantity, the Fourier 
transform of two-level cluster function, known as the Form Factor [2]. It is 
well known that the Fourier transform of the density of states gives important 
information about classical behaviour of the system, i.e. it generates periodic 
orbit length spectrum of the classical system. Similarly the form factor yields 
important information about periodic orbit structure in the classical system. 
Before we elaborate this point, let us derive a semi-classical expression for 
the form factor. We will use definition of |l3j to derive expression for the 
form factor 

00 

62(r) = - 2 y" drr2(r)cos(27rrr) (5.48) 


where r is a dimensionless time, related to t by 

r = — tn = 2'KMav 
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Figure 5.3: Two level cluster function: for complete Tr/S-rhombus billiard 
(see text) 
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Figure 5.4: Two level cluster function: for 7r/3-rhombus billiard pure rhom- 
bus modes(see text) 



115 



Using (|5.47|) the semi-classical expression for the form factor is given by (since 
only finite numbers of levels are considered, we will integrate above equation 
from r = to 600 = Ar/2, this will not change our results since Y2{r) ~ 
for large r) 



&2(r) 



'^f* 1 vrAr 



-7^ 

/a ,max 



[Si((7o 



Svrr) 600) + Si ((7a,max + 27rr) b^)] 



'Ja , max 



sin (27rr6oo) cos (7a,maxfeo 



+ %^si(7„,„.ax&oo) [&00 sin (27rr6oo) + ^^^^^g^] + 

(1 - w) [Si ((7a,min - 27Tt) b^) + Si ((7„,^in + ^TTt) b^)] 



+ ^ sin (27rr6oo) cos (7a,min6) 



^si(7, 



a,mm'-'oo / 



boo sin (27rr6oo) + 



cos (27rTf)oo) 



(5.49) 

where 7„,max(min) = 27r/„,min////,max{min)- As 600 bccomcs large enough, then 
we can approximate sine integrals as follows. 



vr cos(x) cos(x) 
Si(±a;) ~ ±-=F si(x) ~ — 



2 X X 
Substituting above relations in ()5.49|) we get, for r < mina(7Q,^max/27r) 



(5.50) 



1. f \ r'"maxSin(7„,max&oo) - ?^minSin(7a,minfeoo)T sin(27rr6c 



vrAr 



And for r > max„(7„,max/27r) 



(5.51) 



T / \ J '"max COs(7o,max^cxD) '"min COs(7a,min^cxD 

) I cos(27rr6oo) , ^/ 1 x 

02(r) = - 2^ /«a < T \ 7 — )• 

' vrAr vrrOno 



(5.52) 

Various important properties of the form factor has been discussed in |15| 
[Til [T3] . We bring out some new (not reported so far) features of the form 
factor. In general, form factor 62 (t) saturates to for r ^ 1 due to discrete 
nature of the spectrum. For r ^ 1,62(t) —1- In case of chaotic systems 
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diagonal approximation leads to monotonically increasing (~ r) form factor. 
The contribution of the off-diagonal terms however cancels this rise to restore 
correct asymptotic behaviour of the form factor [T3] . 

From ()5.49|) one can see that there will be a spike at r = 0, due to 
sm{27iTb) / irrb (which is ~ 5(t)) factor on RHS. For < r < minQ(7Q,,min/27r), 
62 (t) will be -1 since contribution of the same factor becomes negligible and 
J^af^a = 1- It should be noted that above expressions for 62(t) have inde- 
pendent contributions from the classes of periodic orbits in the diagonal ap- 
proximation. The classes of the periodic orbits are characterized by length of 
shortest periodic orbits (/min,o) that can be arranged in increasing order. For 
r ~ minQ,(7Q,min/2vr), corresponding contribution from that family will de- 
crease(in absolute sense) and saturating ultimately to as r becomes greater 
than 7Q,min/27r, resulting in net increase in the form factor. If lengths of short- 
est periodic orbits in different classes are well separated then r may still be 
less than next (7a,min/27r), the form factor will saturate at — J2a' > — 1, 
where summation is now taken over all classes except one(i.e.a). This thus 
leads to a "step" like structure in the plots of the form factor. If however, 
lengths of the shortest periodic orbits are not well separated steps may merge 
into each other. The step structure is thus a hallmark of existence of multi- 
ple, distinct " characteristic time scales" (hence also length scales) in these 
systems. 

In fig. 5.5 and 5.6 we study effect of variations in the length parameters 
and K^s on the form factor. Again here choice of various parameters is same 
as that in case of the two-point cluster function. Though general conclusions 
that can be drawn are similar to that discussed in the previous section, i.e., 
the existence of larger length parameters or more weightage (via. k param- 
eters) to higher length parameters shows deviation away from Poissonian 
behaviour. Most important feature is imprints of classes of periodic orbits. 
We have considered three classes for all the cases and the form factor rises to 
zero in the same number of " steps" . Existence of steps in the form factor are 
thus signatures of the classes of bands of marginally stable periodic orbits 
in the classical system. This feature is apparent in the plots of the form 
factor in some of the works [T31 IT^ where systems, not far away from the 
integrability, like kicked rotor or domino billiards are considered and in the 
examples that we will consider in this paper. We wish to draw attention of 
readers to the fact that the form factor saturates at in long r limit and 
does not rise monotonically as in case of chaotic system in the diagonal ap- 
proximation used above. This in turns implies that off-diagonal corrections 
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may not be considerable as far as the form factor is concerned. As our study 
indicates the off-diagonal contributions may only affect the form factor in 
the transition region where it rises from —1 to 0. Effect of the off-diagonal 
terms on step like behaviour is difficult to guess. Recent approach(^16j) to 
treat off-diagonal terms in terms of diagonal contribution however, indicates 
that step like behaviour may be preserved even after consideration of off- 
diagonal contributions. Secondly, since diagonal form factor saturates at 
as expected, off-diagonal terms may only change the transition region most 
likely flattening it. 

Our analysis thus brings out a important feature of the form factor that 
has not been paid attention so far. With increasing energy, rise in these steps 
become sharper and sharper taking the form factor closer the Possonian one. 
Also if differences in the length parameters are small, this step structure may 
not be easily visible. This is evident from Fig. 5.7 and Fig. 5.8 where we plot 
form factor for complete rhombus and for even states respectively. Difference 
between two cases is self evident, rise is sharper in the former case than the 
later case. 

In this chapter we have obtained closed analytical expressions for two 
interesting quantities that are used to study spectral fluctuation properties 
using diagonal approximations. This is an important step which should ulti- 
mately lead to complete understanding of spectral fluctuations by including 
other contributions such as off-diagonal terms, diffraction effects etc. 
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Figure 5.5: The form factor: Effect of variation of the length parameters (see 
text) 
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Figure 5.6: The form factor: Effect of variation of the k parameters (see 
text) 
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Figure 5.7: The form factor: For complete 7r/3- rhombus biUiard (see text) 
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Figure 5.8: The form factor: For pure modes of 7r/3- rhombus biUiard (see 
text) 
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lype 


Class 


Closing 
Point 


Band 
Area 






centre-centre 


odd-odd 


(g/2,p/2) 


Ar 


1 

9 












2 
9 


2,/3C 




odd-even 


{(1,P) 


Ar 


1 

36 


v/21£ 






(2g,2p) 


2Ar 


1 

18 


2v/21£ 




even-odd 




Ar 


1 

36 


v/39£ 






{2q,2p) 


2Ar 


1 

18 


2v/39£ 


centre-edge 


odd-odd 


(3g/2,3p/2) 


5Ar 


1 

3 


3^37C 




odd-even 


(3g, 3p) 


^Ar 


1 

12 


3^7C 




even-odd 


(3q, 3p) 


SAr 


1 

12 


3^7C 



Table 5.1: Summary of results for complete 7r/3-rliombus billiard. 
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Type 


Class 


Closing 
romt 


Band 

A 

Area 


Repetitions 






centre-centre 


odd-odd 


(g/2,p/2) 




all 


_8_ 


^J3C 








8Ar 


even 


5 

36 


2^3C 








8Ar 


odd 


1 
4 


V^c 


centre-centre 


odd-even 




^Ar 


all 


2^ 

lo 


^93£ 








SAr 


even 


5 
144 


2^93>C 








SAr 


odd 


1 

16 


V93£ 


centre-centre 


odd-even 




^Ar 


all 


2 
18 


^39£ 








SAr 


even 


5 
144 












odd 


1 


^39C 


centre-edge 


odd-odd 


(3g/2,3p/2) 


12Ar 


all 


1 
3 






odd-even 


(3g, 3p) 


VIAr 


all 


1 

12 


3v/6l£ 




even-odd 


{3q, 3p) 


12Ar 


all 


1 

12 


3^19C 



Table 5.2: Summary of results for 7r/3-rhombus billiard: Only Pure Rhombus 
Modes. Here Ar is area of V, half equilateral triangle. 
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Chapter 6 



Short and Intermediate Range 
Spectral Fluctuation Measures 

6.1 Introduction 

We shall obtain expressions for three important and popular spectral fluc- 
tuation measures namely spacing distribution, number variance and spectral 
rigidity using diagonal approximation for pseudo-integrable billiards. Spac- 
ing distribution provides information about spectral fluctuations on short 
energy scale i.e. on the scale of mean level spacing. Other two measures are 
called as intermediate range spectral measures since they provide information 
about fluctuations on few mean level spacing. In this chapter 

6.2 Nearest Neighbour Spacing Distribution 

Nearest neighbour spacing (NNS) distribution is the probability density P{S) 
of finding the nearest neighbour of a given level at r, in the range r + S to 
r + S + dS. As stated above P{S) distribution of any given system enables 
us to study fine scale texture (i.e., on the scale of mean level spacing) of the 
spectrum. 

For Integrable systems spacing distribution is very well modeled by Pois- 
son distribution {P{S) — exp(— S")) since the energy levels are random due 
to number theoretic properties. On the other hand in case of chaotic sys- 
tems with time reversal symmetry (modeled by GOE of random matrix the- 
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ory(RMT)) spacing distribution is given by the Wigner distribution function. 



provided the level sequence has been normalized to unit mean level spac- 
ing. RMT is expected to be applicable only on those time scales where the 
variables associated with the classical dynamics are random enough, to fully 
randomize the matrices, associated with the corresponding quantum opera- 
tors. Moreover, due to underlying assumptions on which framework of RMT 
is build, it is not expected to shed any light on non-universal behaviour of 
spectral fluctuations which is the characteristic of a given system. 

Various formulas have been proposed to analyse P{S) distribution for the 
systems in the intermediate regime 12 E] ■ However, these models are not 
applicable for the pseudointegrable systems because of basic differences with 
other mixed systems. 

The aim of this section is to develop a model for P{S) distribution in pseu- 
dointegrable systems. Once such model is developed it can be extended to 
more general class of mixed system where almost regular and almost chaotic 
states coexist. When the spectrum is unfolded to unit mean spacing of neigh- 
bouring levels everywhere, one can use probabilistic arguments e.g.|ll|3] to 
develop expression for P{S) distribution. 

The conditional probability g{S)dS of finding a level in the interval (r + 
S,r + S + dS), given one at r is related to P{S) in the following way. Choose 
a segment of length A(integer), and divide it into small intervals, all of the 
same length e. We place A markers at random, independently of one another, 
with the probability g{S)dS, into the small intervals. The first marker above 
will hit any particular small interval with probability £g{^) / A, and miss any 
other small interval with probability 1 — eg{^)/X, where ^ is some coordinate 
inside the small interval in question. Now we choose a contiguous interval of 
length S, and require the probability that none of the markers fall within S, 
while there is a marker in interval {r + S,r + S + dS) . We have to form the 
product of all the 1 — eg{^)/X for < ^ < S, and multiply with eg{S)/\. In 
the limit of small e or large A we find that 



C on RHS is used so as to satisfy condition /q°° P{S) = 1. Underlying assump- 
tions in deriving this relation are: 1) Given a level at r, the probability that 




(6.1) 




(6.2) 
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another level will be around r + is proportional to S and does not depend 
on r. This assumption though applied for all S, is valid only for small S. 
In other words, this assumption means (l)the two point correlation function 
R2{fii^2) is linear in | ri — r2 | and (2) the probabilities in various intervals 
of length S/m obtained by dividing 5* into m equal parts are mutually in- 
dependent. In other words three-point and higher correlations are negligibly 
small. Though both of these the assumptions are inaccurate, above relation 
gives accurate result for GOE of random matrix in terms of Wigner surmise, 
which indicates that errors arising due to these assumption almost cancel 
each other j^]. It is, therefore reasonable to expect that above relation would 
yield good approximation to P{S) distribution in pseudointegrable billiards 
as for such systems levels are expected to be less correlated than that in 
the chaotic systems. Strictly speaking, ^(5*) depends not only on two-point 
correlations but on higher correlations too [HIIZI- 

In case of homogeneous spectrum (i.e. density of states almost inde- 
pendent of energy ), unfolded spectrum with unit mean spacing everywhere 
conditional probability density g{S) is nothing but the two- point correlation 
function, or g[S) =1 — Y2{S) |11IH]- In terms of this the P{S) distribution 
can be written as 



P{S)=Ce~^ {l-Y^iS)) e 



j^Y2(S')dS' 



(6.3) 



In our case C = 1. Using expression for 5^2 (5*) derived in the last chap- 
ter and carrying out necessary algebra one can write NNS distribution for 
pseudointegrable billiards as 



P{S) = A{S) exp{B{S) - S) 
where ^(5) = 1-Y2{S). And 



(6.4) 



B{S) = E«:a{^[Si(7a,max5) + 



sin(7Q 



+ 



COs(7a,max5') 



72 52 

'a, max 



si(7a,max»S') 



vrAr 



Si(7«,min5) + + COs(7«,„.i„5) 

' CK ,min 



+ 



"Si(7a,min»S') 



(6.5) 
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where Si(x) = 7r/2 + si(x) and 7q, 



min(max) 



271"^a,min/^H, 



min(max) 



. The expres- 



sion for ^(5*) is clearly normalized. For small S, neglecting terms of order 
0(5*^) and higher, one obtains 



A{S) ~ As 
B{S) = B 



J2a ^ 



27 ■ 



•\/^max '\/^\nir 



S-AS 



s 



(6.6) 



thus 



P{S) ~ Aexp[-AS] 

Therefore dP{S)/dS \s=o = —A"^ < 0, showing level attraction near S = 0, 
which is weaker than that of a Poissonian distribution for a finite number of 
levels but approaches the same in the semi-classical limit. For large S again 
our expression approaches Poissonian distribution. To see this behaviour 
clearly we study from (j6.4|) . dependence of P{S) on length parameters and 
parameters k s. we will consider only three classes of the periodic orbits. 
Fig. 6.1 shows the effect of variations in the length parameters and Fig. 6. 2 
shows the same for k^s, parameters chosen are same as that in case of 
two-point cluster function. Again conclusions that can be drawn from these 
figures are similar to that stated in the previous chapter. In the high energy 
limit, again P{S) distribution converges to a Poissonian distribution. For 
finite energy levels however, P{S) distribution shows mixed behaviour of 
weaker level attraction as well as level repulsion with respect to Poissonian 
and GOE respectively, for different values of 5*. Note, a closeness of curve 
(d) in Fig. 6.1 to GOE for S < 1, though curve do not rise as much as that 
of GOE one may easily misunderstood it as showing GOE kind of behaviour. 

In Fig. 6.3 and 6.4 we show plots of P{S) distribution for complete rhom- 
bus and even states of the rhombus. P{S) for complete rhombus approaches 
comparatively faster to Poisson distribution reasons are obvious as it is a su- 
perposition of four different modes. In Fig. 6.4 for curve (a) we choose levels 
'"min = 47 and rmax = 370, to compare our results with that of numerical 
results of ^1] for approximately same number of levels. As stated earlier our 
results are in good agreement with numerical results for S < 1. This can be 
more clearly seen from Fig. 6.5 where we show cumulative P{S) distribution. 
The curve (a) in this figure is our results and numerical results are shown 
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Figure 6.1: Nearest neighbour spacing distribution. Efect of variation of 
length parameters (see text) 
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Figure 6.2: Nearest neighbour spacing distribution. Efect of variation of k 
parameters (see text) 
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by triangles and dots (for two different mode of even states). Results are 
in close agreement for S < 1, however for S > 1 deviation of our results is 
too large to be correct. This behaviour of P{S) distribution using diagonal 
approximation is also in contradiction with recent studies([niE]), where one 
get 

P{S) =ASexp{-2S) (6.7) 

Again for S << 1 our results (for finite number of levels) are in agreement 
with eq. ()6.7|) but deviates appreciably for 5* > 1. This contradiction may 
be removed by considering off-diagonal contributions. As the calculation 
of P{S) involves exponential of integration of ¥2(3), any error in Y2 will 
be amplified exponentially and effect of this error will be more for S > 1. 
Further improvements are therefore necessary in (|6.4|) to take account of 
off-diagonal contributions. 

6.3 The Number Variance 

In the preceding section we have discussed one of the important measure 
for spectral fluctuations on the finer energy scales. On the intermediate 
scales or on the scale of few mean level spacings, important measures are 
spectral rigidity and number variance. Both of these statistic provide similar 
information on spectral correlations. In this section we will discuss number 
variance for the pseudointegrable billiards. The number variance S^(-Z^; tq) is 
defined as the variance of the distribution of the number of levels in intervals 
[ro, ro + L] , n(L; tq) = N{ro + L) - N{ro), 

ro)=<[n{L;ro)-L]'>, (6.8) 

where, < • ■ ■ > denotes a usual spectral averaging over tq. A completely 
random spectrum that follows Poissonian behaviour shows a linear trend 
for a number variance, S^(L;ro)=L, whereas on the other hand, for the 
chaotic systems where levels are strongly correlated, a number variance is 
asymptotically given by S^(L; ro) ~ (2/7r^)/n(L) + .44 for L > 1. The number 
variance E^(L;ro) depends on Y2{r) via relation (e.g., see jT2]) 

S2(L;ro) = L-2 {L - s)Y2{s) ds (6.9) 

JO 
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Figure 6.3: Nearest neighbour spacing distribution, complete 7r/3- rhomb us 
biUiard (see text) 
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Figure 6.5: Cummulative nearest neighbour spacing distribution. 7r/3- 
rhombus billiard (only even modes) (see text) 
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Using results of last chapter, we can exactly integrate above equation to get 
simple analytical form for E^(L;ro) as, 



Lsin(7a.,max-L) 



+^^^^oos{^^,^^L) + cos(7„,^axi^) 



6 -""V ia,max^/ < ^-y, 
si(Ta,max-^) 3^ 



TrAr 



+ 



■ COs(7„,mm-t') 



37q 



COs(7a,min-^) + 



L3 



Si(7a,min-^) 



37q 



(6.10) 

where again 7a,min(max) = 27r/„,min/^if,min(max)- For small L, one can see that 
S^(L;ro) ~ L i.e. it shows Poissonian behaviour, as L increases, we can 
approximate Si{x) ~ 7r/2 — cos(a;)/a; and hence si(a;) ~ cos(a;)/a;, S^(L;ro) 
starts deviating from the Poissonian and oscillates around a non- universal 
value as 



TrAr 

2 



Lsm(7c,,maxl/) 

6 



TrAr 



COs(7a,maxJ/) _ 

37a ,max 37a ,max 

Lsin(7a,min^) COs(7a,inin^) 



6 



37 



37q 



Thus the saturation value for very large number of levels (i.e. 7q, ^ L) is 
given by 




Kir 



' max ' min 



Q ^ mill. a 



Ar 



This behaviour is very much akin to that of the spectral rigidity and satura- 
tion value is correctly given by Sqo — 2A300, where Aa^oo is obtain directly 
(i.e. not from Y2) in the sext section. In the Fig. 6.6 and 6.7 we show this 
behaviour as well as effect of variations in the length parameters and k pa- 
rameters respectively as done earlier. Trend is same. Hence our conclusions 
are also similar to the one stated earlier and hardly needs any elaboration. In 



136 



Fig. 6.8 and 6.9 we show more realistic cases namely complete rhombus and 
even states of the rhombus respectively. In Fig. 6.9, one can compare our 
results for energy level windows (47,370) (curve (a)) with numerical results 
of [13j. Discounting oscillations of numerical results agreement is good even 
for r > 1. This is because, number variance involves simple integration over 
I2 and any error due to off-diagonal contributions is not as much amplified 
as that in the P{S) distribution. In the next section we will confirm this in 
the case of another spectral measure known as spectral rigidity. 

6.4 The Spectral Rigidity 

In this section we obtain ^3] the Dyson-Mehta As-statistic for pseudointe- 
grable billiards and show that it is nonuniversal with a universal trend, also 
that this trend is similar to the one for integrable billiards. We present a 
formula, based on exact semiclassical calculations and the proliferation law 
of periodic orbits, which gives rigidity for the entire range of L. To consol- 
idate our theory, we discuss several examples finding complete agreement 
with the numerical results , and also the underlying fundamental reasons for 
the nonuniversality. 

In their statistical theory of energy levels of complex systems, Dyson and 
Mehta [T3] proposed the A-statistic to study spectral fluctuations on the 
intermediate energy scale, the most popular being the As-statistic defined as 
a local average of the mean square deviation of the spectral staircase from the 
best fitting straight line over an energy range corresponding to L mean level 
spacings. 

(/ 7-i\ Lf'^^av 1 
min^^ j de[N[E + e)-A-Bsf^ (6.11) 

— L l'2dav 

where N{E) is spectral staircase introduced in chapter(l) (keeping in conso- 
nance with existing literature we use symbol L for correlation range in this 
section, and not to be confused with symbol for side length used earlier ). 
(...) represents spectral averaging over energy scale much larger than what is 
called as outer energy scale but much smaller than classical scale [16], which 
is already discussed in the previous chapter. Though, A3(L) can also be 
expressed in terms of the two-point correlation function, we prefer to use 
above direct definition to obtain an expression for A3 in pseudo-integrable 



137 




1 2 j, 3 4 5 



Figure 6.6: The number variance: Efect of variation of length parameters, 
(see text) 
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Figure 6.8: The number variance:complete 7r/3-rhombus billiard (see text) 
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Figure 6.9: The number variance: for even modes of 7r/3-rhombus billiard 
(see text) 
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billiards. Minimizing 16. llb ver A and B we get 



ML) 



, , L/2dav 

J deN\E) 



L/2dav 

J deN{E + e) 

L Jldav 



-12 



.2 L/2da 

av{EL\ J 
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edeN{E + e) 



- L l'2dav 



(6.12) 

For a Poisson spectrum (by which integrable systems are very well mod- 
eled) A3(L) = L/15. Consider another trivial example of a spectrum of 
the harmonic oscillator in one dimension. Its spectrum is also called picket 
fence spectrum due to equi-spacing between levels. N[E) for such case is a 
ideal staircase and obviously mean square deviation of a best fit from this 
staircase is constant (1/12) independent of L. For the chaotic systems with 
time reversal symmetry which are modeled by GOE of random matrix the- 
ory rigidity is given by = InL/n"^ — .007. It may be noted that for 
L << 1, the staircase nature of N{E) leads to the limit A3 L/15, what- 
ever distribution(non-singular) the levels have. We now proceed to obtain 
expression for spectral rigidity in psedo-integrable billiards. 

E 

A asymptotic semiclassical approximation for N{E) = J d{E')dE' in case 



of pseudointegrable billiards is given by 



N{E) 



2-^ ,3/2 



cos 



'Ijy/E 

h 



(6.13) 



Substituting this in eq. 16.121 we can obtain the rigidity as 



,3/2,3/2 
k l-j I'k 



G{y,,yk) (6.14) 



where 



G{yj,yk) = fivj - Vk) - f{yj)f{yk) - ^f'iyj)f'iyk) 



(6.15) 

and f{y) = siny/y, y = LI / {Adav{E)E^/'^K) = txLI/Ih- All other symbols 
have their usual meaning. 

The main questions that ensue from the numerical studies [T71 [TTJ [THJ HH] 
are : a) are the levels of PIB uncorrelated and mimic a Poisson process in a 
certain range of L, as seen in IB ?; b) is there any saturation of rigidity in 
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PIB if L exceeds the system dependent range, as seen for the IB ?; c) what 
is the essential difference between PIB, IB and chaotic billiards in terms of 
level correlations ?; d) can we obtain a formula for the rigidity such that 
Poisson and non-Poisson results follow in a natural way for the IB and PIB 
respectively ? . Here we answers all these questions to a large (sometimes 
complete) extent. This success holds due to the fact that the rigidity is a 
direct consequence of the proliferation law with some simple, nontrivial mod- 
ifications in the known formalism shown below. Recalling discussion in the 
last chapter, we employ the uniformity principle ^U] and retain also, apart 
from the diagonal, the off-diagonal part corresponding to the systematic de- 
generacies in the lengths of po's (giving same contribution to the rigidity as 
the diagonal terms). From the exact results on some of the PIBs obtained 
in previous chapters, one can classify the bands of the po's in such a way 
that the projective phase space area occupied by all periodic orbits in a given 
class (defined by corresponding lattice points where po's closes) (say a) is 
identical (Aa). With this in mind, we can write ^^{L,E) as 



A3(^,i^) = SS 

L ^'J 

(6.16) 

where G{y) = 1 — /^(y) — 3(/(y))^ and Greek subscripts again denote classes 
of periodic bands. In ()6.16p . gag 13 denotes the number of po's with the 
same action belonging to the class a(/3), and 6ij is the usual Kronecker 
symbol. In the above equation the summation J2j G{ya,j)/la,j can be written 
as /j^.^ dFa G{ya)/la in the continuum hmit due to the mathematical nature 
of the summand, where dFa represents number of periodic orbits within 
length / and / + dl. This dFa can be deduced from the proliferation law 
(average or asymptotic part) as we have done in the earlier chapters. We 
then have (after unfolding spectrum via rescaled energies tq = Edav) 

AsiL, To) = [Ea 9lAlaah,a + Er, '^i^,/^^7^r, ^ ^S^l,'?] 

+ g^3/2^3/2^i/2 [Ea QaAlpah^a + E7 Er, Sl^,l^g'ygrjAyAnhr,l2,ri] 

(6.17) 
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where 



oo 



2/min,Q: 

and 

oo 

h,a = j dya y~^ G{ya) 

2/min,a; 

both these can be evaluated easily with ymin = 7r-^/-^max- And 

- Ih 



mill. a 



^min.a 

where Ih is Heisenberg length already introduced earlier. For small ymim h = 
27r/15, I2 = 1/9; and for large 7/min, h,a = -f'max,a/7rl/, l2,a = LI^^ J2tt'^L'^. 
For L < min„Lmax,a/7r, 

(6.18) 

For L >> maxo, Lmax,a/7r, 

^3(-^;^o) = 27r3l4R [Ea fi'a^a^a-^max.o + E7 Er; ^l-y,lr,9^9ri^"/^r)0'r)Lmax.,r]] 



+ ^0^7/2^^2^1/2 [J2a 9'a^a^Q-^max,a + E7 J2t, ^l^,ln9j9v^l^r]br,L'^ax,v\- 

Ft 

(6.19) 

It is important to note that minimum and maximum (over a) L^^i^^ cor- 
respond respectively to the longest and shortest (over a) orbits of the set 
containing shortest periodic orbits of different a's. It is the consequence of 
this observation that will lead us to understand the fundamental distinction 
between the spectral correlations of integrable and pseudointegrable billiards. 

Ignoring I2 for the sake of brevity, the formula valid for the entire range 
of L is given by (denoting L^^^^^/nL = lH/T^Ll^\n,a by Aa), 

A3(L,Vo) = ^Ea/^a[A«-iA3-AA5 

-^A«cos(2/A«) - sin(2/A«) + ^A^ cos(2/A„) 
+^fA4 sin(2/A«) + cos(2/A«) - ^^si{2/K)]f^ 

(6.20) 
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where si(x) = — dt sin t and = A'^g'^aa/'n'Aji. Eq. (j6.2(J|) along 
with the hmiting results (|6.18p . (|6.19|) ; which clearly establishes relation 
between A3 and information about the classical periodic orbits such as pro- 
liferation law, band areas, degeneracies in lengths etc. This result applies 
to all integrable and pseudointegrable billiards. To understand the formulae 
better, we now propose to examine some paradigm systems carefully, and 
subsequently compare the results with the known numerical results. 

In this regard, we consider the specific examples of an Incommensurate 
Rectangle Billiard (IRB), the Single Slit Rectangle Billiard (SSRB) and the 
7r/3- rhombus Billiard (RHB). 

Firstly, let us consider the IRB with sides (£,7£), 7 being an irrational 
number. The all periodic orbits fall in a single class occupying projective 
phase space area AApt {Ar = 7^^), except the two shortest periodic orbit 
bands parallel to either pair of sides of IRB. The area of these two shortest 
periodic bands is 2Aji. The proliferation law for the IRB can be easily found 
by employing the ideas of stacking and replication, we get (counting all the 
repetitions of the ppo's), 

Firb{1) = al' + bl = ^ (g 21) 

Using a, A{= 4^4^) in (9) we get complete quantitative agreement with results 
obtained earlier [1^] for L < L^a.^/ir. As observed in [12], the oscillations in 
A3(L) are rather weak beyond the "crossover regime". 

To get the correct saturation values of A3(L), we have to consider 0(/) 
term in eq. ()(j.2H) . Taking account of this, in the region where L < Lmax/^r, 
we get 

A3(^, Vo) = - + ( ^4^1 (6.22) 

The second term is quite small as compared to the first one due to ro~°°^^ 
factor. For L ^ L^^-^j-n^ on the other hand, we have 

which is in very good agreement with the numerical results. 

Our next example is the single slit rectangle billiard, which is a simple 
variation of the barrier billiard already discussed in previous chapters. Re- 
call that this is an example of a PIB whose invariant surface is topologically 
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equivalent to a sphere with two handles (genus, g = 3). The law of prolifer- 
ation is the same as a^l + bj- We can obtain aa,ba for different classes of 
bands in this system in similar manner as in chapter (4). 



With these, for L < minQ,(Lmax,o/7r) = ^Jro/in, 



and for L >> maXa(Lmax,a/7r) = ^Stq/tt, 

'V2 + I3 



A3(L,ro) 



+ 
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,1/2 



(6.25) 



We will discuss these results after we present calculation for yet an- 
other well studied system - the 7r/3-rhombus billiard. This is an almost 
integrable system with an invariant integral surface of genus two. For this 
system tessellation of the plane is not complete and results in more gen- 
eral barrier structure Again, recall that each trajectory from origin 
to a coprime pairs, {q,p) represents ppo ending at c{q,p). We then get for 



L < min„(Lmax,a/7r) = j2ro/37V3: 



L I L?' 

A,{L, ro) ~ + 3,/421/2^3/2 ;T72 (^.26) 



and for L >> maXa(Lniax,o/7r) = \/2-\/37rro/3 




A3(L,ro)~.237-^ + y^ (6.27) 

From expressions and examples discussed above, one can clearly see that 
there is an universal trend of A3(L) with L for integrable and pseudoin- 
tegrable billiards. More precisely, for L < mma{L^s.x,a/'^), the rigidity is 
very well approximated by L/15, and for L » ma.Xa{Lj^sLx,a/ t^) it saturates 
with a crossover connecting these two limits smoothly. The extent of the 
crossover region is given by the difference between miuo, and max^, of L^ax,a, 
or in other words, depends on the spectrum of lengths of shortest periodic 
orbits over a. Nonuniversal aspects, thus, arise due to nontrivial classifica- 
tion depending upon the degree of tessellation of invariant surface in terms 
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of a system-specific fundamental region. For instance, in IRB, tessellation is 
complete and there is only one class of bands (a = 1); the crossover region 
is expected to be of lesser extent - a fact fully corroborated by the numer- 
ical experiments. In the SSRB, there is a barrier (gap to barrier ratio is 
unity) in a rectangle which gives rise to a periodic untessellated arabesque 
in terms of which classification is facilitated - the number of bands here is 
seven. Similarly for the RHB, the number of bands is eighteen. Importantly, 
it should be noted that the value of L at which the spectral rigidity devi- 
ates from the Poisson value of L/lb, and the value at which saturation sets 
in, depends upon the lengths of the shortest periodic orbits distributed over 
various classes admissible in a given system. Indeed, this is the fundamental 
source of nonuniversality. 

Let us discuss the numerical results on various pseudointegrable billiards. 
Most of the studies have been on rhombus billiards JI] , square torus billiard 
and its generalizations 117^ and singular billiards [22] • The analysis for the 
singular billiards was carried out and one understands the level spacing statis- 
tics . The study of the two-level cluster function (in particular T?{L)) 
does not give the GOE result j21] although the level spacing is GOE raising, 
thereby, a question currently beyond explanation. Therefore, we concentrate 
to explain the results for non-singular systems. 

Perhaps the paradigm pseudointegrable billiard is the RHB [TT| ITH] . In 
both these studies, one can observe that the rigidity is intermediate to Poisson 
and GOE. From our analysis, taking the energy and parameters from these 
numerical work, it turns out that deviation from L/15 would occur at L ~ 1 
and 2 respectively. We illustrate this in Fig.6.10, where we compare our 
analytical result with that of the numerical work JT] and the agreement 
is clearly evident. Here also as in the case of number variance errors due 
to neglecting off-diagonal contribution do not affect much. The crossover 
values are also correctly predicted by our analysis. We show behaviour of 
A3 for small values of L in Fig. 6. 11, where deviation from L/15is evident. 
In Fig. 6. 12 we show A3 for complete range of L, where the crossover region 
and the saturation can be seen. Since the numerical results are not available 
for higher energy and for larger range of L, the saturation cannot be clearly 
seen in the numerical experiments. It is, therefore, desirable to carry out 
extensive numerical work for higher energy and for larger range of L. The 
formulae ()6.18l6.19p provide the guidelines for choosing appropriate number 
of levels to bring out all the salient features of the systems discussed above. 
Our analysis also explains the results of T!7] where one gets L/15 for very 
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small values of L and there is a saturation regime. Unfortunately, because 
of constraints over levels available, the belief of an intermediate behaviour 
between that of Poisson and GOE has been pursued for quite sometime. Our 
analysis clearly reveals, that such a behaviour does not exist and the spectral 
rigidity never becomes GOE. 

The occurrence of periodic orbits in the bands is a likely reason for the 
slow rise of A3(L) in large L region and overall stronger fluctuations than 
the GOE result. A recent result on Stadium billiard indicates this possibility 
too ^25J - in this work, /S.^{L) is shown to be rising well above the GOE curve 
if the contribution of the bouncing ball modes is taken into account. In 
chaotic systems like this (also, e.g. the Sinai billiard) the analysis of bands 
can be carried out using the above theory and it is expected that there 
exist a departure from GOE as well as a rise in spectral fluctuations after 
some L decided by the length of the periodic orbits in the band. Recently, 
non-genericity of the rigidity arising from banded orbits is discussed for the 
Stadium biUiard PE] . 

In conclusion, we have obtained good approximation for the Aa-statistic 
in an close analytical form for systems in which periodic orbits of the marginal 
stability (in bands) occur, from which Poisson and non-Poission results fol- 
low in a natural way. Answers to the basic questions we have asked above 
are as follows: a) The levels of PIB are uncorrelated and mimic a Poission 
process for L < minQ,(Ljnax,a) which depends on shortest periodic orbit of a 
given system, hence nonuniversal value. This condition also stipulate mini- 
mum number of energy levels that one should be consider in the numerical 
experiment to observe this effect, b) For L >> maxQ,(Linax,a) which depends 
mainly on the longest of the shortest p.o. among the different classes, the 
spectral rigidity saturates to a nonuniversal value, c) The fluctuation prop- 
erties of PIB and IB differ essentially in the extent of transition region. In 
IB transition region will be of less extent, since there is only one class of 
p.o., deviation from a Poission and saturation is determined by same p.o. (i.e. 
shortest one). In PIB because large number of classes of p.o. are present, all 
shortest p.o. amongst the different classes play important role in determining 
shape and extent of crossover region. 
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Fig. 6.10: Spectral rigidity for /3-rhombus billiard. 
Curve (a) is our result for no. of levels =350, 
triangles represents numerical result and dashed 
line is Possonian L/15. 
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Fig. 6.11: Spectral rigidity for the rhombus. Full 
curves represents rigidity for 350 and 1000 levels as 
marked. \ marked on the curve shows point L=mina 
(L.max, a) where respective curve starts deviating from 

L/15( Poissonian). 
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20. 40. ^ 60. 80. 



Fig. 6.12: Spectral rigidity for the rhombus as in 
fig. 6.11. Start of saturation is marked by \ n 
respective curves at L=max^(Lm^x,c^ 
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Chapter 7 
Conclusions 



The work presented in this thesis provides first step to attain complete un- 
derstanding of spectral fluctuation properties in pseudointegrable systems 
within semiclassical framework. It also supplements our knowledge about 
underlying reasons for universal and non-universal properties of spectral fluc- 
tuations. Furthermore, close analytical forms for various spectral measures 
are obtained, which in our knowledge were not available before this work 
(albeit under diagonal approximation only). 

Semiclassical framework needs complete information about periodic orbits 
of a system under consideration. In particular one needs to specify actions 
of periodic orbits and their stability properties. For example, in the case of 
pseudointegrable billiards we need to know lengths and band areas of the pe- 
riodic orbits in a given system. This is a herculean task, and there are vary 
few systems about which complete information regarding this is available. 
Different methods are developed to gather information about periodic orbits 
in different kinds of system (e.g. symbolic dynamics). These methods are 
some times system specific and sometimes more general. In chapter three we 
have developed a methodology to enumerate and classify periodic orbits of 
some pseudo-integrable billiards. Here we have exploited important geomet- 
rical attributes of systems under considerations. To our knowledge these are 
only pseudointegrable systems about which complete information regarding 
periodic orbits is available. The same methodology can be used to study large 
number of rational polygonal billiards in particular one can apply technique 
to rectangular billiard with many slits(or barriers) or L-shaped billiards or 
its further generalization. This is one of the open problem in which we will 
be interested in future. This will also give us opportunity to deal with higher 
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genus pseudointegrable systems. 

As it turns out from our analysis proliferation of periodic orbits (or growth 
rate) in one of the most important factor in the semiclassical methods. Spec- 
tral fluctuation properties indeed depends on the exact asymptotic prolifer- 
ation of the periodic orbits in the system. We obtained this law in chapter 
four using probabilistic arguments. We proved that the law of prolifera- 
tion of periodic orbits in the pseudointegrable billiards under consideration 
is quadratic in length, in fact , of the form al^ + 6/ + c. We also obtained 
coefficients a, b, c explicitly using probabilistic arguments. 

A complete knowledge about periodic orbits enabled us to obtain ana- 
lytical expressions for various spectral measures. Most important of these is 
two-level cluster function which is studied along with its Fourier transform 
in chapter five. To obtain two level cluster function one have to evaluate 
density-density correlation function first. The semiclassical approximation 
for density correlation function leads us to summation over all periodic or- 
bits. This formulation does not explicitly brings out various factors that 
affect spectral fiuctuations in the system. Using simple but non trivial argu- 
ments we convert this form of summation to the integration. It then turns 
out that proliferation law play important role in this integration and hence, 
in the spectral fluctuations. This conversion to integration enables us to ob- 
tain analytical expressions for various spectral measures. This also enables 
us to study effect of considering finite number of levels on the spectral fiuc- 
tuations. In particular spectral properties depends on l)lengths of shortest 
periodic orbit occurring in different classes of periodic orbits. 2)phase space 
areas occupied by classes of periodic orbits along with respective coefficients 
of proliferation law which forms a factor k in our analysis. In particular we 
show importance of factor lamin/^H, where In = V^tt^TrTo is the Heisenberg 
length. From our analysis it is clear that this ratio for maxo(/Qmin) should be 
<< 1, for one to talk about asymptotic properties {h — > 0) of spectral fluctu- 
ations. For 7r/3-rhombus biUiard this condition indicates that the number of 
levels considered should be at orders of magnitude larger than 1600. For other 
pseudointegrable system this number may still be much higher. Though it is 
meaningful to study statistically large number of levels in the non-asymptotic 
regions, this will bring not bring out universal properties of that class of sys- 
tems as spectral fluctuations are not free of system dependent parameters. 
This chapter thus we are able to deflne classical parameters that governs 
spectral fluctuations and extent of their influence. In the case of form factor 
even diagonal approximation gives correct asymptotic (r >> 1) form unlike 
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that in chaotic systems where it was found that diagonal approximation leads 
to linear rise in the form factor in asymptotic region (r >> 1), which was 
compensated by off-diagonal contribution to give correct saturation. This 
indicates that off-diagonal contribution in pseudointegrable systems can af- 
fect only transition region. Further, we have shown that imprints of different 
classes of periodic orbits can be seen in the behaviour of the form factor. 
This feature is still gone unnoticed though many numerical experiments do 
have these imprints. 

In chapter six we consider short and intermediate range spectral measures, 
namely spacing distribution, number variance and spectral rigidity. We also 
compare our results with numerical experiments available. We found that in- 
termediate range measures are well modeled by our theoretical expressions. 
The departure from numerical results are not serious and may be attributed 
to off-diagonal contributions. Spacing distribution however, agree with nu- 
merical experiments only for short correlation range this is due to fact that 
any error due to off-diagonal contributions is exponentially amplified hence 
rises fast for larger correlation ranges. This can be improved further by 
considering off-diagonal contributions. 

We now enlist some open avenues that our work suggest 

1) To obtain complete understanding off-diagonal contributions should 
be taken into account exphcitly. This is however difficult task. Recently 
developed indirect method by Bogomolny to consider these contributions 
in terms of diagonal contribution may proved to be fruitful in this regard. 
Some preliminary study confirms our conclusions above. It also indicate 
that approach to Poisson like behaviour is further slowed down by these 
contributions. 

2) As mentioned in the thesis trajectories of pseudointegrable billiards 
resides on foliated surface of genus greater that two. This therefore do not 
allow trajectories to explore phase available as in the case of chaotic systems. 
The singular points in the phase space, however, results in diffraction of 
trajectories that hits it. These diffracted orbits are therefore not confined to 
one foliated surface but may visit many more via multiple diffractions thus 
imitating chaotic behaviour. These should be properly accounted for. 

3) Yet, another kind of orbits which are on the edges of the band of 
periodic orbits can give different contribution to periodic orbit sum hence 
should be studied in detail. 

4) Extension of methods developed here to higher genus systems is also 
a one of the important problem. 
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5) Can one generalise these results on two dimensional polygonal billiard 
to n-dimensional polyhedra. 

6) Using semiclassical formalism and complete information about periodic 
orbits can we develop a method to extract information about the eigenfunc- 
tions for these systems. 

7) One of the important applied field where study about polygonal bil- 
liards can be used is semiconductor microstructures where dimensions of a 
system is much smaller that the mean free paths of electrons (hence electron 
imitates billiard like dynamics). To design and develop microstructure with 
desired properties it is important to develop complete semicalssical under- 
standing of such systems. 

I hope this is just a beginning of my persuasion of knowledge about 
dynamical systems in general. 
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